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Kurzfassung 



Das Standardmodell der elektroschwachen und starken Wechselwirkung (SM) be- 
schreibt bis heute fast alle experimentellen Daten innerhalb der experimentellen und theo- 
retischen Unsicherhciten. Zu seiner inneren Konsistenz, dh. im speziellen zur Anomalie- 
freiheit benotigt das SM ein sechstes, zu Beginn dieser Arbeit noch nicht entdecktes Quark, 
das Topquark (t). In den letzten beiden Jahrcn (1994-95) wurde dann die Evidenz fiir 
seine tatsachliche Existenz durch Experimente am Hadronbeschleuniger Tevatron immer 
starker, sodafi heute der Nachweis bei einer Masse von ca. 180 GeV als gesichert gilt. 

Der ideale Platz fiir eine genauere Untersuchung der Eigenschaften des Topquarks ware 
aber ein e'^e^-Beschleuniger. Dort konnte man Top-Antitop Paare {tt) nahe der Schwelle 
erzeugen und so seine Masse, Zerfallsbreite und die Kopplungen zu anderen Elementar- 
teilchen genau bestimmen. Durch die hohe Masse und den raschen Zerfall in und b 
besitzt das Topquark die einzigartige Eigenschaft unter den Quarks, daB sein Verhahen 
rein storungstheoretisch bestimmt werden kann. Die storungstheoretische Behandlung 
der Bindungszustandseffekte nahe der Schwelle fiir die Erzeugung ist Gegenstand der 
vorliegenden Arbeit. 

Eine der einfachsten Erweiterungen des SM ist das sog. Minimale Supersymmetrische 
SM (MSSM). Die darin vorhergesagten skalaren Partner zu den existierenden Fermionen 
konnten in ahnlicher Weise erzeugt werden und werden daher in dieser Arbeit auch 
behandelt. 

Im Kapitel 2 werden zuerst sowohl fiir instabile Fermionen als auch fiir Skalare 
Bindungszustandsgleichungen angegeben, die als Grundlage fiir eine systematische 
Storungstheorie dienen. Als Beispiel wird das Massenspektrum zweier gebundener Skalare 
zu 0{a^) berechnet. Dann wird das tt-Potential in konsistenter Weise zu numerischer 
Ordnung 0{aj) bestimmt. Dabei wurden ein neuer Box-graph berechnet sowie Beitrage 
der elektroschwachen Wechselwirkung beriicksichtigt. Im Rest des Kapitels wird dann 
die Frage einer moglicherweisen grofien Korrektur zum Wirkungsquerschnitt auf Grund 
einer laufendenSZerfallsbreite behandelt. Es kann gezeigt werden, daB ein bisher nicht 
beriicksichtigter Beitrag zu groBen Kompensationen fiihrt. Durch die Einfiihrung der 
neuen Gleichung fiir instabile Fermionen und der Verwendung einer Wardidentitat lafit 
sich schlieBlich die Rechnung soweit vereinfachen, daB der allgemeine Mechanismus dieser 
Kompensationen sichtbar wird. Dies fiihrt zu dem Theorem, daB fiir alle nicht durch 
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Annihilation zerfallenden schwach gebundenen Systeme die Bindungszustandskorrekturen 
zur Zerfallsbreite allein durch die Zeitdilatation beschrieben werden konne. 

Kapitel 3 beschaftigt sich mit der Berechnug der Wirkungsquerschnitte nahe der 
Schwclle mit numerischen Methoden. Zuerst wird der allgemeine Formalismus beschrieben 
und der totale Wirkungsquerschnitt fiir die ti Produktion diskutiert. Dabei zeigt sich 
unter anderem die Wichtigkeit der groBen Zerfallsbreite fiir eine storungstheoretische 
Behandlung einerseits sowie die Bedeutung der Resummierung im Gluonpropagator 
andererseits. Dann wird die Vorwartsriickwertsasymmctrie nahe der Schwclle erstmals 
ohne die Einfiihrung eines kiinstlichen Cut-offs bercchnct. Die dabei verwendete Methode 
lafit sich auch auf die Berechnung des axialen Beitrags zum totalen Wirkungsquerschnitts 
bcim tt-System und die Bestimmung des fiihrenden Beitrags zur Stop-Antistopproduktion 
anwenden. 

Konnten in dieser Arbeit wichtige Aspekte der Produktion schwerer, stark wechsel- 
wirkender Teilchcn nahe der Schwclle geklart werden, so bleiben doch noch cinigc theo- 
retische Herausforderungen fiir die Zukunft. Dabei sei speziell die Anwendung der ri- 
gorosen Methoden aus Kap.2 auf die Berechnung physikalischer Observablen wie in Kap.3 
erwahnt. Dies wird notwendig werden, um Theorie und Experiment mit Genauigkeiten 
von kleiner als 1% zu vergleichen. 
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Chapter 1 
Introduction 



Perturbative expansions in the coupling constant in quantum field theory possess two 
types of applications, the calculation of scattering amplitudes and the computation of 
processes involving weakly bound systems. Many of the successes of quantum electrody- 
namics (QED) are, in fact, related to positronium, i.e. to the second one of the aforemen- 
tioned applications. The proper starting point for any bound-state calculation in quantum 
field theory is an integral equation, comprising an infinite sum of Feynman graphs. The 
Bethe-Salpeter (BS) equation fl]] fulfills this task and it is well known that in the limit 
of binding energies of O(a^m) the Schrodinger equation with static Coulomb attraction 
is obtained. The computation of higher order corrections to the Bohr-levels, however, 
turned out to be far from trivial. It was recognized, though, relatively late that, at least 
conceptually, substantial progress with respect to a systematic treatment results from a 
consistent use of a perturbation theory geared to the original BS equation p|. As an addi- 
tional bonus one avoids nonrelativistic expansions as implied by Hamiltonian approaches 
with successive Fouldy-Wouthuysen transformations 0. Within the BS-technique, how- 
ever, it is desirable to have an exactly solvable zero order equation different from the 
Schrodinger equation, because otherwise e.g. the approximation procedure for the wave 
function lacks sufficient transparency, especially in higher orders and may even break 
down in certain cases, one of which is treated here. One of the advantages of the BS 
approach to perturbation theory is the freedom to select a different zero order equation. 
Of course, in that case, certain corrections included already at the zero level are to be 
properly subtracted out in higher orders. An especially useful zero order equation has 
been proposed some time ago by Barbieri and Remiddi (BR equation Still, one of the 
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most annoying features of all bound-state calculations remains the pivotal rule played by 
the Coulomb gauge. In other gauges, e.g. already the (in QED vanishing) corrections to 
0{a^m) of the Bohr levels imply to take into account an infinite set of Feynman graphs 
0. Only in very special cases, when certain subsets of graphs can be shown to represent 
a gauge-independent correction, another more suitable gauge may be chosen. 
By contrast to QED the vast literature on bound state problems in quantum chromody- 
namics (QCD) adheres to a description of the quark-antiquark system by the Schrodinger 
equation with corrections 'motivated' by QCD [^]. As long as a relatively small number 
of parameters suffices for an adequate phenomenological description of observed quantum 
levels, this approach undoubtedly has an ample practical justification. Furthermore even 
for the relatively heavy fiavors charm and bottom nonperturbative effects can be shown 
to be of the same order of magnitude as the leading perturbative ones. Apart from lattice 
calculations there exists no method to do bound state calculations in a non-perturbative 
way. Only if the non-perturbative contributions are smaller than the leading perturbative 
ones it is known how to include them in a systematic BS pertubation theory |0, If. In the 
case of QCD this requires particles with masses in the range of lOOGeV. Especially the 
recent evidence for the top quark with a mass around 170 GeV seems to indicate that this 
particle is the ideal testing ground for perturbative QCD bound state calculations. The 
cleanest environment to observe the top quark clearly would be a next generation e~^e~ 
linear collider. On such a machine other heavy particles with strong interaction could be 
discovered as well. If supersymmetry is realized in nature it seems very likely that one 
of the scalar partners of either top or bottom could be observed. Since it may well be 
that even in, say, a 500GeV collider those new scalar particles (and also the top quark) 
will be produced near their threshold, a reliable prediction of the production cross section 
slightly above threshold seems very important. The signal of the top quark is expected 
to be best suited for an exact determination of its mass, its width and an independent 
measurement of the strong coupling constant. 

Thus also for this reason a return to more rigorous QCD arguments remains as de- 
sirable as ever. The standard literature on quarkonia (see e.g. 0]) is almost exclusively 
based on nonrelativistic expansions or on the calculation of purely static forces |Ty . 
Moreover, very often potentials with higher order corrections as determined from on-shell 
quarkonia scattering are used In these cases relevant off-shell effects which are typ- 
ical for higher order corrections may even be lost altogether . On the other hand, from 
the point of view of relativistic quantum field theory as elaborated in the abelian case of 
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positronium, a similar, more systematic approach seems desirable, the more so because 
the basic techniques are well developed. In addition, at least in one case, namely the de- 
cay of S-wave quarkonium, the result of a full BS-perturbation calculation including 
the QCD corrections to the bound state wave function, yields a result very different from 
the one which took into account only the corrections to the quark antiquark annihilation 
alone [|l3 |. 

In this context the relatively large size of the running coupling constant even at high 
energies represents a well known problem, together with large coefficients from a perturba- 
tive expansion. Therefore e.g. problems arise in the comparison of the coupling constant 
as determined from scattering experiments within the minimal subtraction scheme (MS), 
with the coupling constant to be used in a consistent weak bound-state approach. The 
philosophy within our present work will be that the orders of magnitude, as determined 
from ajjg will be used for estimates, but that we shall imply a determination of as by some 
physical observable of the quarkonium system itself. In that way delicate correlations of 
'genuine' orders of from basically different types of experiments are avoided. 

From high precision electro-weak experiments of the LEP collaborations, the mass 
range of the top quark now seems to be established to lie in the range lQ0-200GeV 
IM, in agreement with the direct searches at the tevatron |T5[. Thus for the first time 



a nonabelian bound state quarkonium system seems to fulfill the high mass criterion 
required for a genuine field theoretical approach. Unfortunately the drawback of this 
situation is that the weak decay t ^ b + W broadens the energy levels |TB| for increasing 
mass rrit so that above ^ 120GeV individual levels effectively disappear. 

But this disadvantage turns into a virtue since it completely eliminates nonpertur- 
bative and renormalon contributions. The reason for this can be stated as follows. The 
time scale ^qcd necessary to allow strong interactions to become effective and to produce 
e.g. hadrons with open top cannot compete with the decay time Fg ^. Therefore the top 
quark decays before nonperturbative strong interactions can act on it. In this situation 
"pure" quantum chromodynamic (QCD) perturbation theory is sufficient for a complete 
description of a strongly interacting system. This applies as well to Coulombic QCD 
bound-state effects in toponium although also there the hydrogen-like levels are strongly 
smeared out by Fq. Nevertheless, the computation of those levels represents a necessary 
first step since the order of magnitude of these corrections can be used as a order param- 
eter for the determination of the quark-antiquark potential. Furthermore we will see in 
chapter 2 of this work that a lot of insight into the physics of toponium can be gained 
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from the study of the properties of the four point function at the bound state poles. 

After all, this is the first instance where the application of quantum field theory for 
weakly bound nonabelian systems is justified because even nonperturbative effects can be 
taken into account in an equally well-defined perturbative manner. 

A correct formulation of QCD in Coulomb gauge entails not only Faddeev-Popov- 
ghost terms but also the inclusion of nonlocal interaction terms [|^. Therefore, the full 
Lagrangian reads (a=l,...,8 for SU(3)): 

1 ^ 

where the Lagrange multiplier B"^ imposes the Coulomb gauge condition, and where 

= d.-tgT'^A;, (1.2) 
F;, = d^At-d^AI^ + gUcAlAl. (1.3) 



Vi and V2 are given in and will be discussed more explicitly below. The above La- 
grangian will include all effects of the strong interaction, but, as we will show, QED and 
weak corrections may also give contributions within the numerical order of our main in- 
terest. This effects will be discussed within the framework of the Standard Model (SM). 
With respect to possible extensions of the SM we will restrict ourselves to the Minimal 
Supersymmetric SM (MSSM) where we especially focus on the occurrence of new heavy 
scalars like the stop and the lightest (uncharged) higgs boson. 

The thesis is organized as follows. In chapter 2 we first construct a solvable zero order 
equation for decaying particles (sect. (2.1)-(2.2)). Then in sect. (2.4) we calculate within 
the framework of Bethe-Salpeter perturbation theory (sect. 2. 3) the potential for the top- 
antitop [ti) system to numerical order O(a^) . As a further application of these methods 
we derive the bound state corrections to the toponium decay width in sect. (2.5). 

Chapter 3 is devoted to the calculation of cross sections to be measured in a next gen- 
eration linear collider. Within the Green-function approach we study some consequences 
of the results of chapter 2 for the total cross section (sect. 3.1.1). Furthermore in sect. 
3.1.2 we propose a method to circumvent an unphysical singularity in the calculation of 
the forward-backward asymmetry. The same method may also be applied to the calcula- 
tion of the axial contribution to the ti total cross section, as well as to the production of 
heavy scalar particles near threshold (sect. 3.2). 

Acknowledgement: This work has been supported partly by the Austrian Science 
Foundation (FWF), project P10063-PHY within the framework of the EEC- Program 
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Chapter 2 



Application of Rigorous Bound State 
Methods to Toponium 

2.1 A Relativist ic Equation for Decaying Fermions 

In this section we will present a solvable relativistic equation for decaying fermions, similar 
in form to that for stable fermions of Barbieri and Remiddi 0]. While this equation is 
applicable to variety of bound state problems we will focus especially on the tt system. 

The Bethe-Salpeter (BS) approach for weakly bound systems starts from the BS- 
equation 




G = D + D K G 



(2.1) 

for the two fermion Green function G{P,p,p'), where D is the product of the two (full) 
propagators and K represents the 2pi BS-kernel. All four point functions in |2?^ depend on 
the total momentum P = {Pq, 0); the incoming (outgoing) lines carry momentum P/2±p 
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{P/2±p'). 

Near a bound state pole the Green function assumes the form 

Gij,i-j{Pi,P2,Pi,P2) = {2ny6^{P - P)^Xij{p,n,l)- Ti-Xij{p,n,l) + Greg (2.2) 

where XijiP: is called BS- or bound state wave function. 

Perturbation theory starts from an equation similar to eq. with Gq, Dq and Kq 



chosen in such a way that the exact solution is known. If Dq consists of the nonrelativistic 
propagators only and if Kq is the Coulomb kernel the zero order equation corresponding 
to ^]T] (after integrating out po; p'o etc. ) simply reduces to the Schrodinger equation. A 
very convenient zero order equation which already includes the relativistic free fermion 
propagators and still remains solvable has been found by Barbieri and Remiddi (BR) 
some time ago 0]. However if the decay width of the constituents becomes comparable 
to the binding energy of O(a^m), perturbation theory runs into troubles because of the 
occurrence of terms F / (a^m) in graphs like fig. |2.1| , containing a chain of subgraphs which 
are responsible for the decay. To circumvent these difficulties one has to include at least 
a part of the exact self energy function in the zero order equation. 

E E E 

^w,* -^w,* ■^W.D 



Fig. 2.1 

The BS equation in terms of Feynman amplitudes for K and S reads for a bound state 
wave function x 



xrfip;P) = -^Su'{^+p)s,,,{-^+p) J ^ir,,,,,,„(p,p,pOx^?'5"(p';P), (2.3) 



where 5* is the exact fermion propagator, and K is the sum of all two fermion irreducible 
graphs. Furthermore, we have introduced relative momenta p and p', a total momentum 
P = Pi—p2, and we choose the center of mass (CM) frame where P = {Pq, 0) = {2m+E, 0). 
The Green function reads in terms of field operators 

Gijjj{pi,P2, Pi,P2) := J t^'xirf%irf^X2rf%2e*^P^"^+^"'"^-^^"^-^"^"^^Gij;Jj(xi, X2, Xi, X2) 
G,,jj{xi,X2,Xi,x-2) := (0|r^i(xi)^,(x2)^i(xi)%(x2)|0) (2.4) 
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In terms of Feynman amplitudes and the relative coordinates this can be written as 



Gijj]{P,p,p') 



d Xd xd xe 



i(X P+xp—xp) , 



(2.5) 



|Tvl/,(X + -)^,iX - )vI/,(- )vI/.( )|0) 



The notation can be read off the pictorial representation in fig |2.2| . i,j are collective 
indices for spin (cr, p) and color (noted 





Fig. 2.2 

It is well known that the dominant part in K for weak binding {a 0) is the one- 
Coulomb gluon exchange which results in an ordinary Schrodinger equation with static 



Coulomb potential. With the help of the nonrelativistic scaling argument [|18 



Po ~ 0{ma'^), \p\ ~ 0{ma) 
Po ^ 2m-0{ma'^) 

this result is even independent of the chosen gauge: 

Ana 



(2.6) 



K 



-7° ,7° 



(2.7) 



Focusing on QCD we notice the fact that only color singlet states can form bound states 
because the Coulomb potential is repulsive for color octets. The color trace will always 
be understood to be already done, leading to the definition 



^4^ 

" 3 4:71 



4 

-as 
3 



(2i 



to be used in the following, in terms of the usual strong coupling constant as- Because the 
above mentioned nonrelativistic limit of the BR equation contains the projection operators 
A^, defined below, it is awkward to calculate the so-called relativistic corrections in a 



straightforward way within the framework of BS perturbation theory, starting from p]3 
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with Therefore we will use the generalized BR equation as described below instead 
of the Schrodinger equation. 

Moreover, for our present case, we need a generalization of the equation given in 
for unstable fermions, described by complex m ^ m = m — if o/2 where f o/2 represents 
the imaginary part of the self energy graph responsible for the weak decay of the free 
top quark. Since the real part of m should still determine the sum of polarizations in 
the numerator of the propagator, we take by analogy to the stable case in the zero order 
approximation to |2.3| the relativistic free propagator {Ep = i/p^ + m^, rj = mTQ/2Ep) 



5(±Y+Po) 



[(±y+p)7-"^] ^ 
A+70 



+ 



A 7o 



±f + po - + ir] ±f+po-Ep- ir] 



+ 0i 



m 



(2.9) 
(2.10) 



with the relativistic projectors 



j^±^^ ^ Ep ± {dp + pm) 
2E„ 



[2.11) 



If furthermore dKo/dpo = both sides of p.3| may be integrated with respect to po. On 
the r.h.s. the product of the two propagators 5* with p.9| yield four terms, with Cauchy 
poles determined by irj, two of which give no contribution. Generically we obtain 



2txi 



A+7o ® A-70 A-70 ® A+7o 



Po - 2Ep + 2ir] Pq + 2Ep - 2ir] 



so that Ko{p,p')^o with 



^0 = [ Xodpo 
Ztt J 



(2.12) 



(2.13) 



remains to be inserted at the places of Written with indices the direct product 

notation can be explained as: A® B := Aaa'Bpijs. As in |^ the kernel is now chosen as 



i^o^o = [7oA+A+A+]<l>o[A'-A-A-7o]ir 



so as to annihilate the second term in |2.12| . In 

(fp' 



(Po - 2^p + 2iri)% 



i2ny 



K{p,p',Po,iri)A+X+A+%A' A^A" 



:2.14) 



(2.15) 
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instead of A the projector A := A (—p) appears. Thus expanding $o = J2 A'^$q"^'^''A-^ 



A,B=± 



only ■* is found to differ from zero and obeys 



(Po-2i?, + 2^r7)A+$(+-)A- = J ^^k{p,p' , P,,i7^)K+\+<^^^-\p')\-K- . (2.16) 



The nonrelativistic projection of |2.16| onto A"*" (g> A with 



A+A+A+ = — 

= — (2.17) 



2Ep 



E„ + m 



and the introduction of appropriate factors in K relatively to the Coulomb kernel 



k = ^i^i'vv'K^ (2.18) 



,2 



4/(Po + + 2iri) 
and in the wave function 

%{p)(xu-^^i-^(t){p) (2.19) 

also for the case of an unstable fermion lead to a Schrodinger equation for the wave- 
functions ^{p) in momentum space {E = Pq — 2m + 2ir] = E + 2iri): 

The eigenvalues for En = E + ^ in |2.20| clearly occur at the real Bohr levels E^, i.e. 



E„ = M°-2m = -^-^ + 0(a«); = 2mv'l - = _ (2.21) 

In addition to the selection of the " large" components by the choice ^.14| , obviously also the 
sign of ir] in p.l8 was crucial for the dependence of p.20| on the combination E + 2iri alone. 
Still, for the bound-state argument at complex values of the energy, the independence of 



1] with respect to the momentum p is essential. Going back to |2.9| we observe that with 
the choice 

f = ^r, (2.22) 
m 
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where F = const is the (cm.) decay rate of a single top quark, the disturbing p-dependence 
in 7] is cancelled. This leads to 

r 



The full BS-wave function (color singlet, BS-normalized |T9[) for the BR-kernel can 
be obtained by going backwards to Xo^~^ '■= x{p)- K appears as the real residue of the 
complex pole. The wave function {un = Ep — M°/2 — ie) 



Xnip, e) = 70X.(P, -eho = '^r^)^)^,'^^^' 



(2.23) 



is identical to the BR wave-function of stable quarks and belongs to the spectrum of 
bound states P„ = M° — iVt- The ie has been introduced to determine the integration 
around that pole. 

In eqs. |2.23| S" is a constant 4x4 matrix which represents the spin state of the 
particle-antiparticle system: 



S 




singlet 
triplet. 



(2.24) 



is simply the normalized solution of the Schrodinger equation in momentum space, 
depending on the usual quantum numbers (n, /, m) [0, a±i, oq in p.24| describe the triplet 
states. In the following it will often be sufficient to use the nonrelativistic approximations 
of eqs. |2.23| (u;„ ^ u)„ = {p^/m + E„)/2 — ie) 



xivY 



^^mS = -loXiprio 



[2.25) 



A similar calculation is needed in order to obtain the solution for the zero order approxi- 
mation to the BS-equation itself. Written in terms of Feynman amplitudes this equation 
reads: 

zGo = -Do + DoKoGo, (2.26) 
It can be shown that the solution to this equation is 

2uj 



Gn 



i{2n)'6{p-p')Do 



Po-uj^ 



[A+A+(A+)'7o] ® [7oA-A-(A-)'] x (2.27) 



{2n)^6{p-p') Gc{E,p,p') 



2uj 



2J 



l\2 
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where 



E = — , (2.28) 

4m 

CO = E,-^^, (2.29) 



and Gc denotes the Coulomb Green function, given e.g. in ||21|| . 

In the following we will calculate corrections to the position of the bound state poles 
of this Green function. One may therefore ask whether it is possible to write it as a sum 
over all resonances and the continuum. The Coulomb Green function Gc clearly can be 
decomposed in that way, but inserting this spectral representation into p.27| leads to an 
expansion with Pq dependent wave functions 

Go = E^^^^%^^#^ + ^ (2.30) 

n Pq P-n 

with 



i2n)'5{p-p') , [A+A+(A+)-7o] ® [7oA-A-(A-)>^ 



D = i{27T)^S{p-p')Do+ (2.31) 

2uj 

{Po - ^^)^ 
A+SA^ u(v) 

X„(Po,p) = ^^2V2luMp). (2.32) 

(Po - ^ ) ^(P) 

The summation over n represents the sum over all quantum numbers: discrete [n, I, m, S, Ms) 
and continuous {k, I, m, S, Ms). We emphasize that uj now depends on the total momen- 
tum Pq and not on the fixed energy P„. Furthermore also the normalization of the wave 
functions had to be changed due to the different denominator in |2.3CI| compared to ^^2 . 



It is interesting to note that the wave functions |2.32| , in contrast to the wave functions 



2.23| , also fulfill an orthogonality relation for different principal quantum numbers. 

tr[J ^^Xn(Po,P) p2 _°^2 Xm(Po,P)] = S^m (2.33) 

The net result of this subsection is that the width of a particle can also consistently be 
included for a relativistic zero order equation similar to the nonrelativistic case where the 
replacement E ^ E + iT in the Green function has been proposed first [|^. Here we have 
the replacement Pq ^ Pq + iT except in the free propagator for the small components. 
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2.2 An Exactly Solvable Zero Order Equation for 
Scalar-Scalar Bound States 



In the last section we developed the bound state formalism for decaying particles. Until 
today only fermionic matter fields have been discovered, but in supersymmetric theories 
for each fermion two scalar partners are required. Since some of them, probably stop or 
sbottom, could have masses within the reach of the next generation of e~^e~ accelerators, 
even the observation of bound states of those particles seems possible. These objects 
and systems built of scalar composite particles in atomic physics underline the need of 
an equally clear and transparent approach as presented above for the fermionic case. 
A recent attempt in this direction splits the boson propagator in a particle and 
anti-particle propagator to be able to treat them like fermions. The spectrum is then 
obtained by constructing the Hamiltonian via a Foldy - Wouthuysen transformation and 
a perturbation theory a la Salpeter. This approach suffers from several drawbacks. First 
it will break down in higher orders due to the appearance of higher powers in p and 
second the introduction of a "new" propagator raises questions on whether one is really 
calculating consequences of the original field theory. A difference in our result to that of 



ref. p2[ seems to confirm this suspicion. 

To the best of our knowledge there exists no attempt in the literature to construct a 
solvable zero order equation for the BS equation containing two charged scalars interacting 
via a vector field. Thus we will generalize the approach of the foregoing chapter to this 
case. 

2.2.1 Stable Particles 

As starting point we present here an exactly solvable equation for stable scalar particles 
which interact via a vector field and show how to calculate the fine structure of such a 
system within our approach. 

Starting from eq. |2.3| we would like to use beside the relativistic scalar propagators 
the kernel due to the Coulomb interaction 

Ac(p.p') = tea '^°+P°+Po)(^°-P°-rt) . (2.34) 

[p — p'Y 



This kernel has the drawback that it is po dependent and the exact solution of eq. 2.3 



with 2.34 is not known. However, in the nonrelativistic regime by the scaling argument 
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(eq. |2.6| ) we can start with an instantaneous approximation to the kernel since po is of 
O(a^m) in this region and may be included in the corrections afterwards. Then we can 
perform the zero component integration on the propagator 

dpo 1 



—I 



2vr [(f + PoY - El + le] [(-f + p,f - El + le] 



1 



1 



2EpPo 

and it is quite easy to show that 



2Ep - Po 2Ep + Po 



E.iAEl - Pi) 



(2.35) 



Ko{p,p') = 47ra 



(2.36) 



gives a solvable equation with the normalized solutions 



x{p) 

X{P, e) 
to the eigenvalues 



'Ep{Pi - AEl) 



v^[(f + Po? - El + le] [(-f + pof - El + le 
-X*(P, -e) 



Mp) 



Pn = Mi'^^ =2mJl-al. 



(2.37) 
(2.38) 

(2.39) 



Eq. |2.38| is dictated by the requirement that x should acquire the same analytic properties 
as the underlying field correlators 

X{P) = /e^^^(0|r$t(^)$(-|)|P,), (2.40) 

m = /e-^^^(Pn|T<l>(|)<l>t(-|)|P„). (2.41) 

Using the integral representation for the step function which is included in the time 



ordered product, one derives eq. p. 38 



Taking the equation for the Green function ^.2(j| instead of that for the BS wave 



function and using again 2.3€ we find 



-F{p)^^^^^2M1f[p') 



Am 



with 



E 



Am?- 



Am 



(2.42) 
(2.43) 
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and 



Fip) = UP. , I P.'.. r.. , (2-44) 



[(f + Po)2 - + ^e] [(-f + poy - El + 



These solutions can be used for a systematic BS perturbation theory for scalar con- 
stituents, as will be demonstrated in the next section. 



2.2.2 Unstable Particles 



As in the fermionic case it is desirable to construct a bound state equation for decaying 
particles. This can be done by the replacement 



Er, 



El 



iVm. 



:2.45) 



While p.45| leads to expressions for the BS wave functions which contain very unpleasant 
expressions for the particle poles it has the advantage the propagator has the form as 
expected from the phase space of a unstable particle. Furthermore the above calculation 
remains essentially unchanged if we define the square root in 2.45 to be that with the 
positive imaginary part. Only the energy in the resulting Schrodinger equation and thus 



in |2.42| changes to 



The eigenvalues for Pq are 



0,n 



E 



2m\ 



p2 



Am? 



Am 



+ iT. 



m 



n 



(2.46) 



(2.47) 



In the case of the fermions we managed to construct wave functions independent of 
r. This was possible because the small components of the propagator containing Pq — iT 
instead of Pq + iT were projected to zero with the help of an appropriate kernel. This 
cannot be achieved in the scalar case and thus, surprisingly enough, the scalar wave 
functions look more complicated than the fermionic ones. To illustrate this we will present 
here an approach in close analogy to the fermionic case. Let us consider eq. |2.35| . To 
obtain an expression like eq. 2.12| we define the zero order propagator 



DoiPo,p) 
A±(Po,p) 



A+(Po,p)A_(Po,p) 



(2.48) 



2E„ 



±f +Po - E^ + ±f+Po + El - 
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In order to be able to use the wave functions of the Schrodinger equation with a Coulomb 
potential we choose a kernel 

Ko-- 



Ip+P+ ^rn^fpf/ 



with 



p = Ep-i- 



^ 2 

p± = 2Ep ± Po - ^r, 

Us = 2Ep + Po + iV. 
This leads to the eigenvalue equation 

x{p) ^ 



d^v' Ip+p'+ 4?7i \/ pp' 
Epp+p^J (2^V ^sK q^EpEp, 



xip') 



(2.49) 

(2.50) 
(2.51) 
(2.52) 

(2.53) 



which is solved by the (already BS-normalized |T^, |2^) wave function 
zDo(Mf -zr,p)^ 



xip) 



2E„ + M, 



(0) 



^ 2Mi°\2Ep + MD + 2zr(Mr - 2Ep){Ep - if 



r(0) 



(0) 



-Ep X 



x(2Ep - Mf )(2E, + Mf - 2zr)0(p), (2.54) 

X(p) = X(p)- (2.55) 

The latter equation is the same as eq. |2.38| if we use real solutions of the Schrodinger equa- 
tion. It is preferable within the present context because of the nasty complex expressions 
in 



2.3 Perturbation Theory 

In view of the calculation of the ti potential we observe that the BR wave function for 
decaying fermions is the same as for stable ones, which means that, as far as the system- 
atic determination of the QCD potential is concerned, we may just consider perturbation 
theory for "nonabelian positronium" for its determination. Perturbation theory for the 
BS equation starts from the BR equation for the Green function Gbr = Gq of the scat- 
tering of two fermions which is exactly solvable. Dq is the product of two zero order 
propagators, Kq the corresponding kernel. The exact Green function may be represented 
as 

1 oo 

G = Y.Xn'p TTXni' + Greg = Go E(^G'o)^ (2.56) 
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where the corrections are contained in the insertions H. Bound state poles Pn = Mn — iTn 
occur, of course, only for < 2m. It is easy to show that H can be expressed by the 
full kernel K and the full propagators D: 

H = ~K + Ko + tD-^ - iDo\ (2.57) 

Since the corrections to the external propagators contribute only to 0{a^) to the energy 
displacement [|2^ , the perturbation kernel is essentially the negative difference of the exact 



BS-kernel and of the zero order approximation. 

Expanding both sides of equation P3B| in powers of Pq — P„ , the mass shift is obtained 



AM-i^ = (/io)(l + ih)) + {hogiho) + 0{h^). (2.58) 
Here the BS-expectation values are defined as e.g. 

We emphasize the four-dimensional p-integrations which correspond to the generic case, 
rather than the usual three dimensional ones in a completely nonrelativistic expansion. Of 
course, p.59| reduces to an ordinary "expectation value" involving d^p and $(p), whenever 
h does not depend on po and p'o- 



In p.58| hi and gi represent the expansion coefficients of H and Go, respectively, i.e. 

oo 

H = Y.hn{Po-PnT (2.60) 
n=0 
oo 

Go = Y.9n{Po-PnT~^ (2.61) 



n=0 

Similar corrections arise for the wave functions f^: 

X^'^ = {9iho + \{hM''^ (2.62) 

To illustrate the use of this perturbation theory as well as the new zero order equation 
for scalar particles presented in the last section we will present here the calculation of 
the fine structure of two stable scalar particles interacting via a vector particle. Existing 
calculations ||2^ rely on a mix of Fouldy-Wouthuysen transformation and the iterated 



Salpeter perturbation theory. Our present approach is much more transparent and allows 
in principle the inclusion of any higher order effect in a straightforward manner. To be 
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definite let us consider the stop-anti-stop system. We will calculate the spectrum up to 
order a^. 

Since in the zero order equation we have replaced the exact one Coulomb exchange 
2.34| by Kq as given in |2.36| we have now to calculate the contribution of —Kc + Kq to 
the energy levels. This is shown in fig. |2.3|a . With 



{{-Kc)) 



—Ana 



(27r)4 (27r 

-12 

■p ' -^p 



4Pn \ / E„E„r 



(2.63) 



4m 
.2^ 



{{Ko)) 



4m 47ra 



we obtain 



AM, 



((-i^e + i^o)) 



ma 



(2.64) 



(2.65) 



2 ' <f ' 16n4 

The fact that the p- integrations are well behaved and the result is of 0(a;^) proves the 
usefulness of our zero order kernel. 

■ ! • ' 



K 



BR 



D 




transversal gluon 
Coulomb gluon 
transversal or Coulomb 



Fig. 2.3 



The transverse gluon of fig. |2.3[ b gives rise to a kernel 



Aira 



{p + p 



r\2 



if-p'^Y' 



(2.66) 
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Performing the zero component integrations exactly and expanding in terms of the three 
momenta one obtains to leading order (c.f. p5| ) 



AMt = {{Ht)) 



Ana .p' 



9 \ -*9 



-) 



ma 



^10 



8n4 8^3 16n3(/ + i; 



(2.67) 
(2.68) 



The annihilation graph into one gluon vanishes due to the color trace since the bound 
states are color singlets. 

The net result for the spectrum of two scalars bound by an abelian gauge field differs 
from that of ref [Q. The difference can be traced back to the relativistic correction |2.65 . 
We have also checked the derivative OKq/OPq contributing to hi and the X-graphs of fig. 
2^.g for possible contributions. Our estimates only yield contributions to higher order. 

So far we have calculated the spectrum of two scalars interacting with a abelian gauge 



field. It has been pointed out long ago |26| that in the case of a nonabelian gauge field 
further corrections arise due to the gluon splitting vertices. The vertex correction shown 
below in p.4| .d has been calculated in p6| for the fermionic case. 

Here we will give a calculation of the same contribution for scalar constituents. After 
performing the color trace the perturbation kernel for the second graph in p.4|.e reads 



2.4.d,2 



d^k (Pq + po + p'^ -ko){-Po+po+ p'o) 

(27r)4 g*2(^_^2[(£+p+^)2_^2]p 

Performing the /cq integration and using the scaling 



-{pq) + 



{pk){qk) 



Po 2m + 0{a^ 
k 



a^Po 



ak 



to extract the leading contribution in a we can show that 



H. 



2.4.d,2 



g m pq 



Adding the similar contribution from the first graph in |2.4| .d gives 



2A.d 



(2.69) 

(2.70) 
(2.71) 
(2.72) 

(2.73) 
(2.74) 
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This result differs by a factor 4m^ from the fermionic result which is compensated by a 
corresponding difference in the wave functions to give eventually precisely the same result 
as in the fermionic case 

The corrections to the gluon propagator are trivially identical to the fermionic case as can 
be seen by inspecting the only "critical" step, the po integration 2.63| - In view of the fact 
that the result depends only on the angular momentum and not on the spin this seems 
reasonable. Thus the difference in the spectrum of the scalar bound state to 0{a'^) is 
entirely due to the tree graphs discussed above. 
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2.4 The tt-Potential 



In the last section we showed how the BS-wave function for decaying particles may be 
reduced to one for stable particles. Therefore, contributions from the width can only arise 
from graphs with internal fermion lines. We will discuss such a case in detail in sect. |2.5| . 
For the determination of the potential the effect of the width turns out to be negligible. 
Especially we can extract the contributions to the potential from the level shifts calculated 
in for the tree and one loop QCD graphs. Our main philosophy will be the following. 
We estimate the magnitude of an interaction by its effect on the level shift. That is to 
say our potential to O(a^) will contain all terms which would give rise to level shifts to 
that order, keeping in mind that the only level that could be measured is the IS* level. 
With the at present favored top mass value of ~ 175GeV this would be possible if rrit is 
known to a sufficient high precision from other experiments by measuring the rising edge 
of the total cross section e~^e~ — * ti. This will be discussed in sect. |3.1| . The contribution 
of the tree and one loop QCD graphs to the potential to O^a"^) reads 



Vi 




47ra -2 



^[Ei(-rm,et) - ^ + J ln(;^ + et)] (2.76) 



47rr ^ Q 2 rrij 



In eq. p.76| the first two lines arise from the relativistic Coulomb correction and the 
exchange of a transverse gluon. They are, therefore, almost identical to the positronium 
case if one includes the color factor 4/3 in the definition of a (c.f eq. |2.8| ). The only 
difference is the absence of the annihilation contribution due to the tracelessnes of the 
color matrix at the annihilation vertex. The first term in the third line comes from the 
gluonic one loop correction to the Coulomb gluon (fig. ^.4| .a, |2.4| .b). Fermion loops are 
included with their explicit mass dependence in the next term (fig. p.4| .c). The fourth line 
contains the contribution of the nonabelian vertex correction fig. p.4| .d. All other graphs 
shown in fig. P]^.e-h do not contribute to the required order in Coulomb gauge. 
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Fig. 2.4 



In the following section we will calculate contributions to the potential from two loop 
QCD graphs and from the electroweak interaction. These contributions will turn out to 
give effects to the same numerical order of magnitude as the ones above. Therefore, the 
technique often used to calculate a potential to a given number of loops is shown to be 
incorrect. 

2.4.1 Two Loop Vacuum Polarization 

As can be seen from the dependence of the potential |2.76| from the masses of the light 
flavors, the usual renormalization group arguments relying on massless quarks in the 
running coupling constant do not consistently include the effect of 'realistic' quark masses 
in the toponium system, when a systematic BS perturbation is attempted. On the other 
hand, in a full calculation of effects of O(a^) two loops with gluons cannot be neglected. 
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However, in the one loop case finite quark masses yield terms of numerical order O 
therefore it can be expected that in a two loop calculation quark masses will lead to 
corrections of 0{a^). 

Furthermore, since two loop calculations with massive flavors are very cumbersome 
we circumvent these problems, for the time being, by the following argument, which 
also includes the three 'massless' quarks u,d,s. Because of the Ward identity for the 
Coulomb- vertex , it is clear from the theory of the renormalization group that the same 
corrections can be obtained by expanding the running coupling constant with a two loop 
(gluons+u,d,s) input for the latter which provides also the first 'nonleading' logarithmic 
contributions. The beta function to two loops is renormalization scheme independent for 
massless quarks and its two loop part has been calculated some time ago |]30| : 



m = -Pog' - Pig' - 

1 38 

Here n/ is the number of effectively massless flavors and (3{g) is the solution of 

rs dg' 



In 



PigY 



which reads up to two loops 



In 



1 

Wo 



1 _]^^Pi^j\Po + f3ig') 



-g' 



g 



/?o g'{f3o + f3ig') 



(2.77) 
(2.78) 

(2.79) 



(2.80) 



(2.81) 



Considering this as an equation for g = g{q'^) we 'undo' the renormalization group im- 
provement by expanding with 'small' g"^ (x a (cf. eq. |2.8| ): 



33 -2n. 
aU-a— - — ^ln^+ 



+ 



a 



(IGvr) 



:[(33 - 2n^)Mn^ ^ - 9(102 - f n;) In ^ 



(2.82) 



Clearly the one loop term agrees with eq. 2.76| in the limit ^0. That limit, however. 



is not appropriate here because in this way we would loose terms of numerical 0{ma^). 
For the computation of the rest we need the expectation value of (In^ ^)/'f ^- This integral 
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can be done analytically (Appendix A) and the result is: 

1 2 g 2 9 

{-^) = — ^{— + ^2 + / + l + + / + l +7 + ln— f} 2.83 

27m^ 12 am 

with 

2{n - I - 1)\ ""d^^ (2/ + l + fc)! 

With eq. [^.83| we obtain for the mass shift, induced by the leading logs of the two loop 
vacuum polarization of the Coulomb gluon a contribution: 

AM2100P = --^^l27'[^ + ^2{n + l + l) + Sni + {^i{n + l + l)+^ + \n^)'] + 
Ylm^n^ y 12 am 

+ 288(^i(n + / + l)+7 + ln-^)l. (2.84) 

am ) 

In this expression we have set n/ = 3 as suggested by the number of sufficiently light 
quarks. Whether eq. |2.84| really represents the full two loop quark- gluon vacuum po- 
larization, numerically consistent with other terms 0{a^), must still be checked in a cal- 
culation of the Coulomb gluon's self-energy to two loop order in the Coulomb gauge, i.e. 



going beyond the sample calculation in |2g] . We, nevertheless, indicate the corresponding 
potential 

V2100P = -2^y|^ {^^^ " 2n/)2[^ + 2(7 + ln^r)2] + 9(102 - yn;)(7 + ln/ir)|(2.85) 
2.4.2 QCD 2-Loop Box Graphs 

It would be incorrect to extrapolate from the QED case the absence of corrections to 
0{a^), other than the abelian tree graphs because gluon splitting allows new types of 
graphs. Therefore we will discuss in this section possible sources of new corrections. 

Our first example of a QCD box graph is fig. p.5| .a. Between the nonrelativistic 
projectors of the wave functions p.25| the perturbation kernel from this graph can be 



written effectively as 

dH d^k — to + m){p2 — ko — m) 



- iH2.5.a = I2ig 



(27r)4 (27r)4 [(pi - t)2 - m^][{p2 - k)^ - m^] 

X ^ —P{q,k,t) (2.86) 

{t-kfk\q-k)H\q-^^ 
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with 



P{q,k,^:= {q~k)k 



[{t- k){q - k)][k{t - k)] 



(2.87) 



(t - kf 

After performing the integrations over the zero components t° and k^ we can use the 



scahng argument p.6| to extract the leading contribution. This is justified a posteriori by 
the (infrared) finiteness of the remaining terms. The resulting expression will thus only 
depend on q: 

dh cPk 1 



2.5.a 



-P{q,k,i) (2.88) 



(277)3 (277)3 fc2(^_ ^)2f2(^_ ^-^2(f _ 

The trick to simplify this expression is to write the inner products in terms of quadratic 
expressions that cancel partly the denominator [|. 

1 



P{q,k,^ 



r2 



■[2q'-k'-iq-ky-t 
t\q- kf , P{q- t)2 ^ fc2(g - kf 



k^{q-i)^' 



+ {k-r)' 



(2.89) 



(fc-t)2 (A;-t)2 (fc-t)2 (A;-t)2 

While the original integral |2.88| is finite, the integrals one obtains from a single summand 
in 2.89 are not. But they are easily evaluated with the help of known results within 
dimensional regularization A convenient check of the lengthy computations is that 
all divergent terms must cancel in the final answer which can be written as 



H. 



2. 5. a 



8177 
'128"' 



'12 



77 



(2.90) 



The nonabelian box graph fig. 2^. a gives therefore rise to an O(a^) correction which leads 
to a slight enhancement of the attractive Coulomb force. This is the only contribution of 
this type in our renormalization scheme which is defined as follows: is the pole mass, 
Zciuon and Zip are defined by a subtraction at fi. However, the usual predictions for 
scattering measurements are in the MS scheme. In a calculation of the potential in this 
scheme a host of contributions of the form of eq. p.90| are to be expected. Since they are 
not known yet, and the numerical coefficient in the result p.90| is small, this at present is 
more of theoretical than of phenomenological interest. However, it provides an example 
of a nontrivial contribution to 0{a^) of a two loop graph which certainly is not included 
e.g. in the running coupling constant. The result p.90 is new, but a qualitative estimate 
was mentioned already in [^, ||10| and p6 |. 

^ I am grateful to Antonio Vairo for bringing this trick to my attention 
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-1 





Fig. 2.5 

Considering the graph |2.5| .b we perform the same steps as for the graph 5. a. One notes 
that either the spatial gamma matrices at the vertices or the zero component integrations 
give rise to two additional powers in the spatial momenta and consequently in powers of 
a. Therefore, this graph contributes only to higher order. It is, however, a spin dependent 
interaction. 

Graphs like fig. |2.5| .a and |2.5| .b with crossed Coulomb lines (fig. |2]^.d) are irrelevant 
because they vanish due to group theoretical factors. 

By connecting the quark lines between the interactions in graphs |2.4|.d and |2.5|.a as 



shown in fig. |2.6| one obtains graphs which could be named nonabelian Lamb shift graphs. 
We have checked whether these graphs could give rise to further 0{a^) effects - but we 
did not find any. Thus all these graphs will contribute at the usual Lamb shift level 
O(a^lna). 
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Fig. 2.6 

The 'X' graph in fig. |2.5|.c can be checked more easily for possible new contributions. 



As in the calculation of fig. |2.5| .a it can be simplified to give 

o- 6 /■ dH d^k (fco+P? + m)(tO+pO-m) 

— itL'iKr = Oiq I —Z TT ~ ~ — 



1 , [Aq-t)] 



2' 



The integration over k yields a divergence l/\q\ if q and p tend to zero. Contributions 
within the order of interest may only result from possible poles after the t-integration. 



For simplicity we consider the part of eq. p.91| from the factor one in the second line: 
It = 



dH +pl + m 



{2ixY[{t + p2f -w?]t''{q-ty 
—2 /"i"^' yq^ — xp2 + p^ + m 



dx dy 
Jo Jo 



(47r)2 Jo Jo {yq — xpY — yq"^ — x{p'^ — m?) 



im /"i , f^^^ . X 



dx / dy — — — + 0{a) 



{AnY Jo Jo x'^m? + y{l — X — y)q 

In M + 0(a). (2.92) 



Thus the part of graph f^.5| .c, specified above, has a leading term from 

^''^ 32(47r)2|g1m m 

as g* — > and therefore contributes to 0(a^lna;). The second part of graph [275| .c gives a 
similar result with a different numerical factor. 



We conclude that - opposite to the QED case |3T|]- box graphs may be important to 
0(q;^). They could, in fact, be responsible for changes of the zero order coupling of the 
Coulomb gluon. As a matter of fact, indeed substantially different values of that coupling 
seem to be required in phenomenological fits of lighter quarkonia. 

As illustrated by the explicit calculations above, to 0{a^) beside abelian QED type 
corrections [Bl, RTI], a host of further non-abelian contributions can be foreseen. 



2.4.3 QED Correction 

As a rule, the consideration of electromagnetic effects in QCD calculations is not necessary, 
but at high energies the strong coupling decreases, and in the case of toponium we expect 
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to be of the same order as aqEo- We may obtain this contribution by simply solving 
the BR equation for the sum of a QED and a QCD Coulomb exchange. This results in 
the energy levels 



Po = MO = 2m^l - i^^±^^^ ^ (2.93) 

^ 2m - m4 - "^""t"^' - - ""H""^' + 0(-% 

An^ 2^2 64n4 An^ ^ v 

where Q is the electric charge of the heavy quark, i.e. 2/3 for toponium. Clearly even the 
'leading' third term can only be separated from the effect of the second one to the extent 
that a{ii) and aqEoilj) can be studied separately with sufficient precision. 

2.4.4 Weak Corrections 

While also weak interactions usually can be neglected in QCD calculations, this is not 



true in the high energy region, because the weak coupling scales like \/Gpvn? ^ becoming 
comparable to the strong coupling if the fermion mass m is large. Even bound states 



through Higgs exchange are conceivable [R^. Therefore, we have to consider weak correc- 



tions and especially the exchange of a single Higgs or Z particle, assuming for simplicity 
the standard model with minimal Higgs sector. 
The Higgs boson gives rise to the kernel 

- lErnggs = -i^Gprn^ „ ^ „ ^ i\f2G pvr? —^—^ , (2.94) 

— mjj q'^ + mjj 

in an obvious notation. 

Since we do not know the ratio am/rriH, which would allow some approximations if 
that ratio is small, we calculate explicitly the level shifts by transforming into coordinate 
space. As in Appendix A, we express the Laguerre polynomials in terms of differentiations 
of the generating function, do the integration and perform the differentiation afterwards 
to obtain 

AMHiggs = -m^-j=—Ini{an), (2.95) 

valid for arbitrary levels and Higgs boson masses, with 

am 
nmH 
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and 



InliO'n) 



2n 



E 

fc=0 



n + l + k\(n — I — 1 



k 



k 



[2.96) 



As an illustration some explicit results for the lowest levels are given in tab. 4. 



n 


1 


2 


2 


3 


3 


3 


I 








1 





1 


2 


Inl{o.n) 




a2(2+a2) 




a|(3+6a2+a|) 






(l+ai)2 


4(1+02)4 


4(l+a2)4 


9{l+a3)^ 


9(l+a3)e 


9(l+a3)6 



Tab. 4 



It is evident that eq. p.95| will give a contribution of order Gpm'^ct^ if the Higgs mass is 
comparable to the mass of the heavy quark and should therefore be taken into account in 
a consistent treatment of heavy quarkonia spectra and the related potential to numerical 
order O(a^). We will thus consider also a term corresponding to p.94| 



'^Higgs = -V2GFm^ 



-mjjr 



(2.97) 



correction in our potential. 
Next we consider the contributions of the neutral current, the single Z-exchange and 
Z- annihilation. For the Z-exchange in the SM we have 



Hz = -V2GFml[Yivf-af^5)] 



rni. 



Yi'^ivj - a/75)] 



(2.98) 



with 



Ti -2Qfsm^ e^, 



-3 ) 



(2.99) 
(2.100) 



where T/ is the eigenvalue of the diagonal SU(2) generator for the fermion /. If / is the 
top quark then T/ = 1/2. Because we cannot expect am to be much smaller that mz we 
use the exact expectation value of the Yukawa potential p.96| . may be dropped in the 
Z-propagator since this provides at most a further correction of numerical 0{a): 
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AM, = -j2G,ml(^ J /"" (2-101) 

2v27r o nm^ Zm^n-^ 2 nm^ 

S*^ = 5'(5' + l) is total spin of the quark- ant iquark system. Therefore this expression gives 
rise to a singlet-triplet splitting within the order of interest. 

Z- annihilation may also yield a sizeable effect. The corresponding energy shift is easily 
evaluated 

AMs=o = ^ 3 Sio 2.102 

2 2 

AMs=i = ^ AMs=o-lO-'AMs=o (2.103) 

a J- — m| 

and also produces a singlet-triplet splitting. It should be noted that the two last contribu- 
tions yield singlet-triplet splittings which are as important as the usual Breit interaction 
(cf. the first two lines in p.76| ). The corresponding contribution to V is 

Vz = v^G^m|a5{^[l-^-(5^-3^)(^ + ^)-(5^-^)] + 

^ir),^ 11 



2.4.5 Schwinger- Christ-Lee Terms 

As mentioned in chapter 1, nonlocal interactions have to be added to the Lagrangian in 
Coulomb gauge. We are not aware of any previous attempt to look whether these terms 
give contributions to bound state problems or to some effective potential. 



By analogy to the second ref. |21] we calculate the vi term to 0((7^ 

. 9 
16 



Vl = -g^— I dhd^r'd^r''Al{r')K,j{f-f')Kjk{r-r'')Al{f'') (2.105) 



KM ^ 



4vr|p| 

This corrects the gluon propagator by 



(2.106) 
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In momentum space 6G can be calculated by using dimensional regularization to give 

^GKg,?') = (27r)^.5'^^5(g-gO^^^(-'^™. + ^), (2.107) 



which means that we have a mass shift with the same structure as the one transverse 
gluon exchange (cf. |^), further suppressed by two more orders in a. 

Since the second term f 2 also represents a correction to the propagator of the transverse 
gluon, it can be estimated by the same method to contribute only in higher orders of a 
as well. We thus find that both terms can be neglected even including terms of 0{a^). 

Finally we will collect the results obtained in this section. Both, the large mass and 
the large width of the top quark provide a new field for rigorous QCD perturbation theory: 
In contrast to the lighter quarkonia, (unstable) toponium is a weakly bound system, to 
be treated by Bethe-Salpeter methods in a systematic manner. The large width even 



further reduces the effects of confinement. As shown first by Fadin and Khoze F can 
be included in the (weakly bound) Green function at the threshold in a straightforward 
manner by analytic continuation of the total energy into the complex plane. 

We first showed that a similar trick may be also applied to a different zero order equa- 
tion, the BR-equation. On the basis of this result we describe how to obtain the proper 
potential V for such a Green function, rigorously derived from QCD in a perturbative 
sense. Although a fully 4-dimensional formalism is used, which especially also allows 
the inclusion of off-shell effects, our result allows an interpretation as a correction to the 
Coulomb potential: 

V = V1 + V2 (2.108) 



with Vi from eq. 2.76| and 



V2 = - — (3--)--^^|27^[-^ + 2(7 + lnH^]+576(7 + lnH_ 
9 r 47rr •' 3 



(2.109) 



+v/2G.m,a,^[l-^-(5 -^—^)(^, + ^)-iS -—^)] 

It has been obtained by checking contributions up to numerical order 0{a'^) to real energy 
shifts, calculated independently from the weak decay t ^ h + W . The abelian relativistic 
corrections and the one loop contributions are collected in Vi. 
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The first term of V2 is tlie nonabelian box contribution, followed by the two loop 
corrections to the coulomb gluon propagator. Here the effects of the quark masses are not 
known yet. On the other hand, we show that certain electroweak effects (QED,Higgs,Z) 
are numerically important to 0{a'^). Their respective contributions are listed in the last 
two hues of |2l09 . 



Within a rigorous field-theoretical philosophy it would be incorrect to add, say, a linear 
term to p.76| in order to describe confinement. At best p. 76 could be supplemented by a 
piece oc {G^)r^ which mimics the tail of confinement effects by gluon condensate |]^, |36| . 

In the derivation of our potential we have not only used the level shifts, but also have 
described in much detail new closed forms for such shifts etc. The reason for that has been 
that on the one hand we hope to have given new useful methods to be applicable also for 
future treatments of nonabelian 0(Q;^)-effects. On the other hand certain computations 
of level shifts may be useful in conjunction with semi-phenomeno logical approaches to the 
lighter quarkonia. However, even in the case of the bottom quark the nonperturbative 
corrections become quite large and their perturbative inclusion seems rather dangerous. 



38 



2.5 Relativistic and Gauge Independent OfF-Shell 
Corrections to the Toponium Decay Width 

As already mentioned in the previous chapters, the top quark decays almost exclusively 
into W~^b. Clearly the tree level SM decay width has to be corrected by QCD and weak 
corrections. QCD corrections to 0{as) ||3^ and the leading electroweak corrections 
have already been calculated. The width is also strongly influenced by supersymmetric 
corrections |3^. For alle these aforementioned corrections it suffices to consider the t- 



quark as a free on-shell particle since it is off-shell near threshold only to O(a^). This 
means also that near threshold not only the 0(q;^) corrections to the free decay are needed 
but also bound state corrections to the decay width play a role. The calculation of these 
corrections will be the subject of this chapter. Concerning the O(a^) corrections to the 
free decay only a partial result is known yet |^ . 



2.5.1 Narrow Width Approximation 

In this section we will calculate off-shell and relativistic bound-state corrections for the 
decay t b + W to 0(ag) making full use of the Bethe-Salpeter formalism for weakly 



bound systems As a first step, we will use the original BR equation as zero order 
equation. Therefore the calculations in this section are performed in the narrow width 
approximation. However, we will convince ourselves in the next sections by using the BR 
equation for decaying particles developed in sect. 2.1 that this actually gives the correct 
result. 

We are able to take into account all terms to that order in a systematic and straightfor- 
ward manner. One of the previously not considered contributions is obtained by explicit 
calculation and cancels precisely gauge dependent terms which appeared in an earlier 



45| off-shell calculation. We will see that important cancellations also determine the 
gauge-independent part. 

As proposed first in ref. the bound state effects for the inclusive reaction ae+e-^tt 



near the tt threshold could be taken into account by including the width Fq for a freely 
decaying top quark |^ in the absorptive part of the vacuum polarization, through the 



replacement of the total energy E by E + zTq. Subsequently this approach was applied 
in other detailed numerical studies involving potentials with the running strong coupling 
constant for the short distance Coulomb term and also including a confining part [|^, 
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45| , However, as noted recently in the introduction of a constant decay rate Tq 
in the Breit-Wigner resonance part of this formahsm leads to problems with unitarity. 
Therefore a "running" width r(p^) for a virtual top was proposed with p"^ ^ mf. As 
expected, a naive off-shell extrapolation of Tq leads to ambiguous and gauge dependent 
results, r(p^) was determined indirectly from the vacuum polarization function making 
certain approximations . Furthermore it was claimed that the running width would 
give rise to effects of the order of 10% to the total cross section. Qualitatively the main 



physical effects from t decaying in toponium seem to be quite well understood by now |47 



On the one hand, F decreases by time dilatation and by reduction of phase space, caused 
by the decay below the mass shell. Cancellations should occur with an increase of F due 
to the Coulomb-enhancement induced by the Coulomb interaction of the relativistic b- 
quark, because this overlap effect is known to be important from the qualitatively similar 
muonium system [^. The purpose of this chapter is to show that existing quantum 
field theoretical technology, well tested e.g. in the (abelian) positronium case, may be 
used also here to solve the problem of a gauge-independent "running width", including 
the qualitative effects listed above. Even the Coulomb enhancement can be taken into 
account properly in this context, which to the best of our knowledge has not been possible 
before. 

The correction to the decay width is obtained from the imaginary part of the expec- 



tation values, receiving contributions from the graphs in fig.p^. 
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- m) 



t b p, = P/2 + p 

Pj = - P/2 + p 
-» • 

t - i = m) 



a: Weak off-shell coirections { E ) 



b i 



b ! 



b: Gjulomb enhancement ( A ) 

Fief. 2.7 



It is important to note that eq. |2.58| is not a power series in as- In p.58| it will 
be essential that O(a^) corrections arise from all three terms denoted here if the BR 
wave-functions are used in lowest order. It should be emphasized that starting from the 
Schrodinger equation instead, the inclusion of the proper relativistic corrections to the 
free fermion propagators would have necessitated the inclusion of complicated higher order 
terms in h. We also anticipate that in the present application of this formalism for the 
imaginary part of p.58| , —T/2, renormalization effects need not be considered. 



The graph fig. 2. 7. a gives rise to the BS perturbation kernel 



i[S(pi) (g) {^2 -m) + {i>i -m)® S(p2)](27r)^(5(p - p') 



(2.110) 



where Pi = ^ + p and —p2 = ^ — p' denote the four momentum of the quark and the 
antiquark, respectively. The direct product refers to the t ®t spinor space. The factors 
(jf^i — m) and {^2 — fn) compensate the superfluous propagator on the line without S. 

Since we are calculating higher order effects we neglect the mass of the bottom quark 
although it may be included in principle. For the electroweak theory we use the i?^-gauge 
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with the gauge fixing Lagrangian ( M denotes the mass of the W-Boson) 



M 



,+ |2 



because it ehminates mixed W — cj) propagators. The gauge parameter ^ will not be 



fixed in the following. We obtain for the imaginary parts of S in |2.110| (s = smOw-, 
P± = (l±75)/2): 



— ImS 



where (for > M^, > M'^/p^) 



(2.111) 



m 



2M2 

(p2 _ ^2)2 



(2.112) 



2p4 ' '-^^ 2^ ' 2p2 ^ 

In evaluating the expectation value according to the rule |2.59| we encounter for both terms 
the same trace (cf |2.11| and p.23| ) 



T 



-tr[7oA-StA+7oP±j^^iA+5A-(j^^2 - m)] = 
hp, + uj)[m{l - 4) + (Po - ^)] + Oiat). 



2 ' 

It is seen to gives rise to (spin independent) relativistic corrections of the same order as 
the off-shell ones from E. Inserting p.ll2| into p.lll| we may expand in powers of {p^ — m?), 
because terms of order {p^ — m^)^ already yield negligible orders in a^: 



- ImS(^) 



C 



327rs2 



2 + —)tiw?) + C- — 



(2 + 





2t(m2)— + 2p(m2,0 



+ 0((p'-m2)^) (2.114) 



Here ImE*^'^^ denotes ImS with j) m and without V±. Thus the total contribution of 
fig.p?7[a to the decay width becomes (/iq = -ff^|po=A/„ „) 



-2Im((/io^)) = -2((1 - - -a^)ImS) 



To - 2rocr„ + 



IGtts^ 



4m2 4 



(2.115) 
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where Fq is the well known zero order result which is twice the decay width of a free 
top quark A factor 2 arises here and in the following, counting both self energy 

contributions of the t and i line: 



-2ImS(^)| 



m] 



m , 
2M2- 



[m 



M 



2\2 



(2.116) 



The expectation value ( ^ ~r ) , left over after the integration, is taken in the ordinary 
sense, i.e. between the (normalized) Coulomb wave functions (f)nim{k) in |2.23| . 

In the graphs shown in fig. |2^.b to the required order it proves sufficient to evaluate 



the leading contribution by applying the nonrelativistic wave functions |2.25| . After a 
straightforward but tedious calculation we obtain for the imaginary part of the vertex 
function A^^^ (^ > M^m'^ and 1 > MVm^): 



ImA 



(W) 



327rs 



(2.117) 



m 



M 



1 + - + 3-^) + 4p(m^ i). 



2M2 



In the BS expectation value corresponding to the first term on the r.h.s of p.58| the 
fco integration to this order allows the simple replacement —>■ —uj. With the spinor 



trace (—3) after including the nonrelativistic projectors in |2.25| we arrive with ImAg at 
{q = p-p') 



IGvrs^ 



I At: a 



(2.118) 



Remembering that p = 0{as), the last expectation value of eq. |2.115| can be rewritten 



2 — / ~ 2/ = ~ 



, 2 Ana 
m q"^ 



(2.119) 



where the Schrodinger equation with reduced mass m/2 has been used in the last equality. 
We note already here that with eq. |2.119| the gauge dependent terms proportional p{rn? 
exactly cancel in |2.115| plus |2.118| ! 



So far only the first term in |2.58| has been considered. The second order terms in h 
for the leading correction O(afa) can be summarized as 



r4 



-2Im[((/.«%i/.^)) + {{h^.g.h^^''))] 
-2Im[((/^«^^))((/.f)) + ((/.o^))((/.?^^))] 



(2.120) 
(2.121) 
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For gi (cf. eq. p.61| ) in p.l20| we find it convenient to write 

X ® X 



9i ■-- 



Gq — 



(2.122) 



Po=M„ n 



We first evaluate the contribution r3^'' of the Green function Gq expanded around 
Po — Mn = 2m Vl — cr'^ — 2m^l — ~ m(cr„ — o"^). For our purpose it suffices to calculate 
the leading nonrelativistic contribution since relativistic corrections are suppressed by 
an additional factor O(a^). Therefore we can use the nonrelativistic Green function 
Gnr(yk", k', a), i.e. essentially Schwinger's solution Gnr{k", k', a) to the Coulomb problem 



pT| , supplemented by appropriate factors (feg — cu^) ^ so as to include also the zero 



component of relative momentum. For Pq 7^ ^n.,o thus a ^ an the first term of eq. 
2.122| yields a contribution (cf. p. 251) 



r^^) = -aJ dkdk'dk"xnr{k"){kl + u';;){ihnS) Gnr{k", k' , a) h^^^{k', k)xnr{k) (2.123) 

where again — (/cq + ^n) represents the inverse nonrelativistic propagator in fig. |2]^.a. 
Carrying out the kQ integrations leads to {dk := d^k/{2TiY etc. ) 



^(1) 



2ro / dkdk'dk"^*{k")[l 



m(cr 



]Gnr{k",k',a)h^''''{k',k)<l>{k), (2.124) 



4tu" 



if we use the fact that for Hq'"^ to O(a^) the QCD corrections in H from Kc — Kbr, 
Coulomb gluon exchange etc. are real and independent of k^. 

The factor proportional to 1 from the square bracket in |2.124j exactly reproduces 
the second part of p.l22| , to be subtracted out before the limit a ^ cr„ is taken. The 
limit ((T^ — (T^) just selects the corresponding pole term 0*0„[y(cr^ — cr^)]^^ in 
Gnr{k", k', a). Thus from ^120| only 



To 



r3 = ^(/^o^^^)(-) 

Z UJ 



(2.125) 



remains, {h'^'^^) could be taken from refs. p8| , p6| . Its explicit form will turn out not to 
be relevant here. 

For the first term of 2.121 /if, i.e. the derivative of 2.11CI| with respect to Pq at 
Pq = Mnfl is needed. Again the factor {h^^^) provides additional powers of and hence 
allows the use of the mass shell condition Pq — > 2m in 



= (27r)425(p-p')[7oS®(^7o) + (^7o)®7oS]. 



(2.126) 
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Thus from this term 



f) = -2((/.^^-))Im((/.f))^-^ 



UJ 



(2.127) 



follows. 

In the second part of eq. p.l21| /if'"'^, i.e. the first derivative of H'^^^ , will only 
receive contributions from H = —Kq + Kbr to O(a^) . The instantaneous Coulomb 
kernel Kc does not depend on Pq, nor does the vacuum polarization ||2^ correction which 
had provided the leading contribution to {h^^^) above: 



^(2) 



(2.128) 



Inserting the explicit (real) expression for the BR kernel Kbr-, (eq. |2.14| with F 
obtain 

-10(^)1 



we 



(2vr)= 



2Pn 



2 



and thus 



.(2) 



(2.129) 



(2.130) 



Before adding all contributions we note that precisely the contribution Fi in p.ll5| 
has been quoted already in ref. as a result as a specific off-shell extrapolation of the 
free decay width. It was obtained there by the intuitive replacement ^ — > y/p^ in the 
on-shell expression. However adding F2 of p.ll8| , with p.ll9| all gauge dependent terms 
are found to cancel. Moreover a striking feature of the sum Fi + F2 is that even the 
gauge-independent terms precisely sum to zero as well, which is by no means "natural" 
in view of the explicit form of F and C, as compared to Fq of p.ll6| occuring in the second 
term of p.ll5| . In a similar manner also F3 + F4^'' cancel each other (cf. eq. |2.125| and 
2.1271), leaving only F4^^ in the final result 



boundstate 



Fnfl 



(2.131) 



Recalling that 1 — cr^/2 ^ (1 — (J? /im?)Y^'^ , the residual effect may be interpreted as 
an approximation of the 7-factor for the time dilatation for the weakly bound top. The 
origin of that term, the relativistically generalized Coulomb-kernel Kbr gives some credit 
to this interpretation. For the complete cancellation of the other terms we will present a 
convincing explanation in the next section. 
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Among the cancelled terms the other expected effects can be more or less clearly 
isolated: Fi (fig. |2]^.a) includes the off-shell corrections and the decrease due to the phase 
space reduction. But this is here completely cancelled by the "Coulomb enhancement" 
(fig. pr7[b). It should be emphasized, though, that bound-state effects of the bi or bb 
systems are, of course, irrelevant for total energies near the ti threshold. Our systematic 
perturbative approach also correctly describes the toponium decay to 0{aweak) as a sum 
of the decays widths t —>■ b + and t — > 6 + W~ (cf. the sum of the terms in fig. |2.7|) . 

A reader used to phenomenological calculations may wonder why no running coupling 
in Us or even a more "realistic" potential involving a confining piece has been used. The 
simple answer is that any generalization of this type could in principle completely ruin 
the strictly perturbative approach (and also the built-in gauge-independence) advocated 
her, where e.g. the gluon vacuum polarization (cf. h^^^ above ) must be accounted for 
as a separate perturbation and must under no circumstances be mixed with higher order 
leading logs in perturbation theory, as summarized in a running ag- However, we will 
show in the next section, that this limitation can be circumvented in the case at hand. 

Anyhow, since the real part of the weak bound-state corrections can be estimated to 
be small, the net effect of our result |2.131| for the discrete levels of toponium allows a 
simple interpretation: It consists of the (level- independent) free decay width Fq and a 
very small imaginary correction factor 1 -|- iFo/2m to the Rydberg energy. 
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2.5.2 Including the Width in the Zero Order Equation 

In this chapter the bound state problem for a fermion-antifermion system is considered 
including a finite decay width of the constituents already in the zero order equation as 
outlined in chapter 2. We focus especially on the ti system for which we reconsider 
the explicit calculation of the bound state corrections to the toponium width of the last 
section, which was performed in the narrow width approximation and needed the use of 
second order Bethe-Salpeter perturbation theory. We find that one obtains the same result 
already in first order BS perturbation theory if one uses our present approach. This also 
explicitly shows the dependence of any result at a fixed order of BS perturbation theory 
on the chosen zero order equation. 

Furthermore the extensive cancellations of gauge dependent terms turn out to be a 
consequence of a Ward identity. This cancellation mechanism is shown to be valid for 
general fermion-antifermion systems in the next section. 

The correction to the decay width is obtained from the imaginary part of the expecta- 



tion values, receiving now contributions from the graphs in fig. |2.8| subtracted by the zero 
order propagator and from fig. ^.7|.b. 



m 



m 




p, - p/2 + p 



p, - - P/J + p 




■'W,<D 



Weak off-shell corrections ( D"^ ) 

Fig. 2.8 

Corrections to the fermion propagator have to be included via the perturbation kernel 



(1) 



(iD 



-1 



(2.132) 



Po=M„-ir 

We emphasize that the perturbation kernel now has to be evaluated at Pq = Mn — iT. 



With eq. |33 and /^o"'|p„=M„-.r — Dq Ip^o we find 



hi^^ = + ^^7o) ® {^2 + (A - m) ® (^(p^) - ^ 7o)](27r)^(5(p - p') (2.133) 
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where pi = ^ + p and —p2 = ^ — P denote the four momentum of the quark and the 
antiquark, respectively. The direct product refers to the t t spinor space. The factors 
{-^ — m) and {j)' — m) compensate the superfluous propagator on the line without S. We 
have neglected terms that will contribute to 0(a^r^/m). 

Now we observe that the terms in |2.133| which include E have already been calculated 



in the last section except for the imaginary parts occurring in the trace 
T := itr[7oA-5U+7oP±^A+5A-(j^^2-m)] = 



^(po + uj)[m{l - + (Po - ^) - (2.134) 



and in 



p — m ^ 2mpo — p — ma — irriT. (2.135) 



We will be able to neglect these imaginary parts in the following since they either give 
rise to real corrections ( which we are not interested in) or to imaginary ones from ReS, 
both of which are of O^al^^^jS). Of course, these corrections would be required a in 
complete calculation of the toponium decay width to this order in order to obtain a 
gauge independent result. This goes beyond the scope of the present work. Furthermore, 
these corrections are state independent, expected to be very small and therefore of minor 



practical interest, anyhow. As can be seen from |2.134j the axial part of S does not 
contribute within the expectation value. 

Therefore, the contribution from first order BS perturbation theory to the decay width 
to O(a^r) reads 

((4'))) = ri + r3 (2.136) 

with Fi as defined before 

Fi = -Alm{{-iJ:{pi)®{:^2-m){2n)^5{p-p')))= (2.137) 
= To - 2Foa^ + ^[(2M^ + m')^^^^^ + 2p(m^ 
and a new contribution from the inverse zero order propagator 



,F „ a' 



F3 = -Alm{{-0{^2-m){27r)'6{p-p'))) = -To{l + f) (2.138) 

Because of the inclusion of F in the zero order equation now in the sum F1 + F3 the leading 
Fo terms cancel and the resulting correction is already of 0(q;|Fo). 
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Now we turn to the calculation of the vertex correction (fig. |2.7| .b). This correction 
remains unchanged. However, in sect. |2.5.1| we calculated ImAo := limg^o,p2^m2 g) 
directly. Here we will show that it is possible to use a Ward identity to determine this 
correction. This also provides a good check for our explicit calculation. 

We observe that 

KiP.Q)l_^, = -QQCDT'^^np) (2.139) 

holds at least at the one-loop level. This identity is believed to hold to all orders if one 
takes into account the equivalence of Coulomb gauge calculations in the usual and in the 
background field method and the Ward identities derived in the latter approach [Q. 

In the graphs shown in fig. |2]^.b to the required order it proves sufficient to evaluate 
the leading contribution by putting the fermions on the mass shell and applying the 
nonrelativistic wave functions |2.25|. Using the explicit expression |2.111| for ImS we obtain 



with the help of eq. |2.139 



7oF (2.140) 



327rs2 

P = 1 + :tt7^ 1-— 1 + 3— +4p"^^0• 
This result agrees with that of p^.ll7| and one obtains a correction to the decay width: 



The fact that the gauge dependent terms in the sum Fi + r2 cancel is thus traced 
back to the identity |2.139| . The gauge independent contribution from first order BS- 
perturbation theory is now different from that of |^ (actually it gives now the net result) 
in accordance with the different zero order equation, used in our present context. 

To complete our calculation we finally check the contributions of second order pertur- 
bation theory. It can be shown that due to the fact that the first weak correction now is 
already of order ©(a^F), this corrections do not contribute to the required order. Thus 
our present approach considerably simplifies the calculation of bound state corrections to 
the decay width. 

Summing up all contributions we get the result 

^boundstate = Tl + 1^2 + Ts = ""'^0^ (2.142) 
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in agreement with our previous calculation in the last section and 



|4T| 



2.5.3 A General Theorem on Bound State Corrections 

The above considerations suggest that the result |2.142| is only a special case of a more 
fundamental statement, with broader range of applicability, which we will now derive . 



Consider a fermion whose decay can be described by the imaginary part of a self energy 
function which will have the general form (in any covariant gauge in the relevant sector 
of the theory) 

S(p) = +/^Sv.(p2) +75Sp(p') +;^75Sa(/). (2.143) 

Within the expectation value p.59| the pseudoscalar and axial vector parts vanish. Fur- 
thermore we will drop all factors 0(p^ ~ Z^^)) which is a valid approximation in all cases 
where a^m <^ m — /i. If this is not the case, either the whole calculation makes no sense 
because then one would have to include the entire rung of Coulomb interactions of the 
remaining particle with the heavy decay products or the decay width becomes very small, 
or both. Let us therefore consider the correction to the decay width which results from 
the remaining parts of the self energy in an on-shell renormalization scheme ( m is the 
pole mass of the fermion). We expand S around that mass shell: 

S(p) = z(ImS5 +/^ImSy) + (S'^ + j^S'v,)(p2 - m^). (2.144) 



Terms of 0(Q;^g„^) are understood to be neglected (cf. remarks after eq. p.l35| ) and we 
denote 



Ex := Sx(m2), X = S,V 
d_ 

From eqs. |2.139| and |2.144| we can easily calculate the vertex correction, say, for the upper 



particle line: 

AoA+ = ^T"7o[amSy + 2m{T.'s + mS'y)]A+ (2.145) 
After some algebra we can write to the required order 

hf^ = (ImSs + mImEy + -70) ® {^2 - m) - iD^^^q [ilmEy + 2m{T.'s + mll'y)] 

(2.146) 

h'f^ = — ^[ilmSy + 2m(S'5 + mS'y)]7o ® 7o 



50 



Since 



r = —2ulmEu 



_^ = -ImEsim^) - mlmSy(m^) (2.147) 



and («-Dor=o ~ ^o,r=o)x = we obtain for the bound state correction to the decay with 
to O(a^a^eafc) 

ATtoundstate = -4ir((A- ® {^2 " m))) = -^0^ (2.148) 

Now it is clear from the preceding argument that it apphes not only to the toponium 
system but also to all other systems which can be described as weakly bound fermion- 
antifermion system with unstable components. Another example which is known for a 
long time is provided by the bound muon system where the above result has also been 
obtained first by explicit calculation This result is easily reproduced from |2.148| by 
replacing Fq by the free muon decay width and m by the muon mass. The effect of the 
reduced mass (in our case m/2, negligible for the bound muon ) is entirely included in 
the expectation value. We conclude that the formalism developed here and especially 
the use of the identity |2.139| simplifies the problem of the bound state correction to 
the decay width in a profound way. It is now possible to clearly isolate the underlying 
cancellation mechanism which automatically gives a gauge independent result which can 
be interpreted as time dilatation alone. However, because we may now also include the 
scalar part of the self energy function and since Kq can be any kernel, the result |2.148| is 



found to comprise |^ and |^ as special cases of a more general theorem: The leading 
bound state corrections for weakly bound systems of unstable fermions (with decays like 
t —>■ b + W~^, fi~ —>■ e~ + Ue + i^fi) are always of the form |2.148| . Among the consequences 
we especially note that p.l48| may be safely applied to toponium even if one uses e.g. a 
renormalization group improved potential. While we have presented a field theoretically 
consistent picture of the "smearing" of the individual levels in toponium below threshold, 
it must be stressed that the treatment of the entire threshold region cross section also 
requires an extension of our study to the "continuum" Pq > 2mt for which the existing 
BS formalism must be adapted. 
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Chapter 3 

Heavy Quark Production in e^e 
Colliders 



In the last chapters we calculated the potential and bound state corrections to the decay 
width with rigorous field theoretical methods. These methods will be necessary to compare 
theory and high precision measurements. However, for the near future it will probably 
be sufficient to apply existing less rigorous methods ||3^. In the case of quickly decaying 
particles near threshold the Green function method ^ ^ seems appropriate. Within 
this approach the four point function is approximated by a nonrelativistic Green function 
fulfilling an ordinary Schrodinger equation with an appropriate potential. This method 
works to 0{as) since the respective corrections to the decay width and potential as well 



as the hard corrections decouple to this order |5T|, At O(a^) this decoupling can be 



expected to break down and thus it turns out that the calculation of the cross section 
of decaying particles to order of 1% requires a deeper understanding of the relativistic 
effects near threshold. A completely satisfying and practically tractable method does not 
yet exist. 

In this chapter we will investigate the production cross section of a tt pair near thresh- 
old. Whenever possible we will use the results of the foregoing chapters. 

The cross section for the process e+e~ — >• X, where X is a set of final states X„ can 
be written as 



52 



with 



(3.2) 



If we neglect initial state interaction and calculate the matrix element to the first order 
in the electromagnetic coupling we can write for the s-channel: 



M= (X|T|e+e-) 



(3.3) 



2sc 



9' 



P^-Ml + iMzTz 



(X|jf)(0)|0) 



The SM values for the vector and axial coupling of the Z-boson read 



(3.4) 
(3.5) 



where Tj? is eigenvalue of the third component of the weak isospin operator. For the top 
quark we have Tj^ = 1/2 and Qt = 2/3. 

We are only interested in the contributions of one flavor to the total cross section. 
Investigating the top quark contribution, we have to replace the currents which in general 
include a sum over all kinds of particles by the current involving only the top quark field 
operator \l/(a;). Splitting the currents into vector and axial vector parts we get 



X) 



iqteJu{x) 



2sc 



2sc 



J^{x) = : ^(x)7^*(x) :, 4{x) = : ^(x)7^75^(x) : 
where from now on we abbreviate 



s 
c 



sin 6u, 
cos 9,, 



(3.6) 



(3.7) 



Substituting eq. p.6| into eq. 3.3 we get for the matrix element M 



M = (27r)^5(P-t-t')t^„(t')(-^e2)^[(a(^) + &(^^5)(X|J^(0)|0) (3.8) 



+ (a(^)+6(^)75)(X|J^(0)|0)lMi) 
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where 



(V) VtVeP^ (v) VtaeP , . 

^(^) ^ OjtvP^^ ^(^) ^ _ataeP^^ ^^^^^^ 
z = - Ml + iMzTz, (3.11) 



and t, t' denote the momenta of electrons and positrons, respectively. Since the mass of 
the top quark is much larger than the mass of the Z-boson we can neglect the width of 
the Z-boson in |3.11| . This means that we can treat the a^^^ and 6^^-* as real valued. 

For the production of heavy scalars the above formulas can also be used , but only 
the vector current is present and has to be replaced by 

J^{x) = I : (<l>*(x)9^$(x) - d^^*{x)^{x)) : (3.12) 

where $ represents the scalar field for a given mass eigenstate. To be definite we will 
identify this scalar field in the following with the lightest scalar partner of the top quark 
as predicted by the minimal supersymmetric SM. For a given value 9t of the L/R mixing 



angle, the vertices for the lighter particle can be written as ieQ{pi + P2) with 1^3 



= -Qt, (3.13) 
= {cos^et-2qtsm^ew)/i2sc). (3.14) 

The coefficients a^^^ and b^^^ have thus to be replaced by 

tts = -qt + — , bs = , (3.15) 

Zscz 2scz 

respectively. The lepton trace for positrons with polarization p' and electrons with polar- 
ization p [rrie —>■ 0) can be readily calculated to yield 

L^" = i5p[(l-pW7^(«^''^+&^''H5)(l+P75)7V(«^''^ + &^''H5)te (3.16) 
= [(a(^)a(^) + 6(^)6(^))(1 - pp') + (a^^'W^ + 6Wa(^))(p - p'W^f + t^^ - t'^Ugn ' 
+ 6(^)a(^))(l - pp') + (a(^)a(^) + 6W6(^))(p - p')]e'^^^Xtp 

This simplifies in the CM frame to a purely space like tensor. Let us now focus on the 
total cross section. Since P* is zero in the CM frame, the hadronic tensor Tj^ has to be 
proportional to Sik and therefore the antisymmetric part in L vanishes in the calculation 
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of the total cross section. However, this part of the lepton trace will lead to a forward 
backward asymmetry as discussed below. It is reasonable to expect that a future linear 
e'^e~ collider will work with polarized electrons and unpolarized positrons. Therefore 
we will assume this scenario in the following. The generalization to the general case is 
obvious from |3.16| . 

For unpolarized positrons (p' = 0) the symmetrical contraction of L simplifies to 

= [aWa(^) + 6W6(^) + (aW6(^) + 6Wa(^')p]p2 (3.17) 

The modulus squared of the hadronic matrix elements can be simplified according 

to 



j:i2nr5iP-Px){0\j^mX){X\j^i0m = 2ahsT^f (3.18) 



X 



where 



(27r)4 (27r) 

9,u)G'''{P) (3.19) 



IjiPp ^ \riXY 



and the indices X, Y refer to the vector and axial currents as defined in B.6| and B.12|. Here 



we made explicit the connection to the four point function G{P,p,p') defined in p.5| . Due 
to parity conservation the mixed V-A currents vanish to all orders in the strong coupling 
constant. The electroweak interaction may well introduce parity non-conserving effects 
but near threshold the heavy particles should be almost on shell and thus such effects 
should be suppressed at least by one order in the electroweak coupling for the total cross 
section. Thus we can split off the total cross section into a purely vector and axial- vector 
contribution. 

a = av + aA (3.20) 

For fermions the vector part is the dominant one near threshold, since it allows the 
production of S- waves (with spin S = 1). The axial contribution is suppressed by two 
powers in a^. 

For scalars we have only the vector contribution, but in the nonrelativistic limit the 
coupling proportional to spatial components (in the CM frame) selects I = 1 states. 
Thus the production of scalars is also suppressed by two powers in a^- 
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Collecting all ingredients from above and multiplying by a color factor = 3 we have 
for the contribution to the total cross section from the flavor under consideration : 



with 



72n 

= -js^[cx + dxp]lmG''''{P)a^-^ (3.21) 
727r 

as = -j^[cs + dsp]lmGs{P)a^f, (3.22) 



c^:=(aW)2 + (6W)2 = + 6^ 
dx := 2a(^)6(^) 4 = 2a,6, 



and 



47ra|^mi) 
3P2 



^MA = (3-24) 



3.1 Top - Antitop Production Cross Sections Near 
Threshold 

3.1.1 Total Cross Section 

For the calculation of the total cross section for ti production in the threshold region the 
large width of the toponium system plays a crucial role. It has the disadvantage that 
individual levels disappear since the width becomes comparable to the Is — 2s splitting. 
While corrections to each pole, which lie now in the unphysical sheet, as described above, 
are in principle possible, such a calculation becomes untractable in practice. However, 
as already mentioned in the introduction, the large width provides an effective cut-off 
for nonperturbative corrections and renormalons. This will be demonstrated explicitly in 
the following. To circumvent the calculation of corrections separately for each level and 
afterwards summing those to obtain the Green function, it is convenient to calculate it 
numerically and use eq. p.21| for the total cross section. 



As a first approximation in the region near threshold it should be possible to replace 
the exact Green function by the nonrelativistic one. This Green function fulfills the 
Schrodinger equation 

{-— + V(r) -E~ iT}G(r, r') = -5(r- r') (3.25) 
m 
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For the potential V{r) we would like to use the potential |2.108| since it is e.g. a sys- 
tematic expansion in the strong coupling constant and thus is expected to lead to gauge 
independent results. However, the inclusion of higher order corrections will lead to some 
problems, as we will see below. 



The leading contributions to p.20| arises from the vector coupling (e.g. X = V in 



3.21|) since it produces tt pairs with angular momentum zero (S-waves). Therefore we are 



interested in Z = solutions of eq. 3.25. Due to the symmetry of the Green function and 



the requirement of regularity at r — > cxd and at the origin, the general solution of eq. |3.25 
can be written as 

oo I 

G{r,r') = j:9i{r,r') J] Y:^in')YUn) (3.26) 

1=0 m=-l 

^ J -^>^' J -^<^' (3.27) 



<' > 



• ly^ fy I ly^ • ry^ ry 



where g<{r^) and g>{r^) are regular solutions of the homogeneous equation 

d"^ 1(1 + 1) 

{— - miE + tT- V{r)) - ±^}g{r) = (3.28) 

at r = or r ^ oo, respectively. The actual behavior of the solution g<{r^) (and also of 
the irregular one) for r — > depends on the potential and on /. Due to the presence of 
the width iT the behavior of gy{r) for r — > oo is given by 

Jiim ^>(r) oc e""-'' (3.29) 

with 

a_ := ^m^~E + ^/W+f^. (3.30) 

This exponential damping behavior is the origin of the infrared cut off mentioned already 
in p5i. Consider for simplicity the case E = 0. Then we can write 



where Eb = mo? jl and are the Bohr energy and the Bohr radius of the system, 
respectively . This formula clarifies that if the width F is approximately equal to Eb the 
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form of the potential is only "tested" up to the order of magnitude of the Bohr radius. 
The situation is improved for E < —T since here we have 



a:'<J^, E<-T. 
V m\E\ 

This relation also holds also for small F and thus the Bohr radius again becomes the 
relevant scale. 



The matching condition for and is provided by the 5 distribution in eq. |3.25 



-m = g^{r)g'^{r) - g'^{r)g^{r) (3.32) 

Numerically it is possible to obtain the regular solution at the origin directly by imposing 
suitable boundary conditions. For the singular solution we use the following method 
10, Suppose you have two solutions: 5'<(r) as above and U2{r) an solution determined 
by an arbitrary boundary condition. Then the solution g^ is given by 

g>{r) = c[u2{r)+Bg^{r)] (3.33) 

with 

B = -\im (3.34) 
(7<(r) 

The constant c can be determined from eq. |3.32| . For / = and a Coulomb-like potential 
at the origin we arrive with the boundary conditions 

g^iO) = (3.35) 
9<iO) = 1 

and arbitrary ones for U2 at the result 

G«(f.O) = J!L!!M±Ii^, (3.36) 

Airw r 

w = U25'< - 9<U2. (3.37) 

According to eq. |3.21| the total cross section is proportional to ImG(0,0). With this 
knowledge it is possible to calculate the total cross section near threshold for any given 
potential V{r). We have written a numerical routine for "MATHEMATICAL^" which 
has been checked to give the correct answer for the known Green function of a purely 
Coulombic potential. 
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Before presenting some numerical studies, a short remark on the couphng constant to 
be used seems in order. It is customary to take the strong couphng constant as in the MS 
scheme. EspeciaUy experimental determinations are always given in terms of a'^^^{Mz)- 
But in bound state calculations it is natural to use an defined differently as we did in the 
foregoing. Clearly it should be possible to relate the two schemes. To avoid complications 
from heavy fermion masses we treat 4 flavors as massless. The difference in the bottom 
quark contribution can be expected to be smaller than the experimental uncertainties in 
Us- Therefore we ignore it in the following consideration. The quark- ant iquark potential 
in the MS scheme can be extracted e.g. from to 0{as)- 



Ana 
q ^ 



^ , g2 31 lOn/, 
a(37r/?oln^ + — -— ^, 



(3.38) 



Comparison with p.76| for n / light flavors gives 



«"^ = «M^(l-^(31-^)) (3.39) 

where BS denotes our bound state scheme. Numerically this means that our a is slightly 
smaller than usual which is advantageous for our perturbative calculation. For = 
180, a^^^{Mz) = -117 ± 0.05 get for a renormalization at = l/rB(/i) 

/i = 18.5±lGel^ (3.40) 
a^ifx) = 0.22 ±0.01 (3.41) 
a(/i) = 0.19 ±0.01 (3.42) 

Remember that the above values for a, a^^*^ differ by a factor 4/3 from as- 

For the comparison of the two schemes the terms of the form a^/q"^ were important. 
To compare to relative 0{a1) terms up to a^/g^ are needed in the potential. While the 
only term of this form in our scheme has been calculated in our present work (eq. |2.90| ), 
we are not aware of a calculation of the terms of the order constal/q'^ in the MS scheme. 
Therefore, it is at present possible to compare the strong coupling constant in these two 
schemes only to 0{as). 

Let us first compare fixed order perturbation theory with the renormalization group 
improved one. This corresponds to an investigation of the effect of the gluon loops given 
in |2.76| and p.85| . The leading term is of the form (7 ± In /ir) . This term will appear in 



any higher order. Therefore we add to the pure Coulomb exchange successively terms of 
the form (7 ± In /ir)" coming from the leading log of the n-gluon loop. 
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The results are shown in fig. p.l] . The uppermost hne corresponds to a pure Coulomb 
potential while the line slightly below has been calculated including the contribution of 
the bottom quark. For all other curves the length of the segments is a measure for the 
number of logarithmic terms (e.g. from top to bottom: n = 2,4, oo, 5, 3, 1). 

R 




Fig. 3.1 



Thus we conclude that the inclusion of the effects of the gluon loops is important. There- 
fore the strict application of perturbation theory is not very useful. Instead one is forced 
to sum all leading logs from the gluon loop corrections to the Coulomb line. Here one 
encounters the problem that the resummed gluon propagator has a nonintegrable pole 
and its Fourier transformation is thus not well defined. 

We will impose the following recipe to obtain a resummed potential [|. We calculate 
the corrections perturbatively in momentum space, perform the Fourier transformation 
for each term and then sum up all terms proportional to (7 + log fir). This leads to a 
geometrical series which serves as the definition for the potential for values of r where the 
sum does not converge. 

The pole of the geometric series lies at 

1 



'pole 



, (ll-2/3n^)Qs ■ 



/i 



(3.43) 



^Another possibility frequently used in the literature is to use the two loop renormalization group 
improved potential 
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We expect any regularisation which removes this apparently unphysical pole to influence 
the potential at distances of that order of magnitude. For a top quark of about 180 GeV 
we flnd 

Tpole ~ (3.44) 

which is about two orders of magnitude larger than the inverse damping constant for this 
mass (eq. |3.31|) if the renormalization scale is chosen as 

fi = — (3.45) 

rB 

or smaller, which is also suggested by the momentum scale relevant for the bound state 
problem. Therefore, we can safely sum the gluon loops and we will see that the determi- 
nation of the cross section is only sensitive to quantities well below the pole. 

Since the resulting potential is also the solution of the one-loop renormalization group 
equation for massless loop quarks it is also one loop renormalization scale independent. 
The contributions from the nonleading log's in p.85| and from the mass of the c, b quarks 
are small and can thus be treated in a perturbative manner. Since we have taken only 
the leading logs from a deflnite set of diagrams it is clear which contribution has already 
been taken into account and which one not. 

Let us flrst briefly compare the numerical approach with the standard bound state 
approach where individual bound states are considered and corrections are calculated in 
a systematic way to the position of the pole and to the wave functions as described in 
sect. p73| . 
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Fig. 3.2 



In fig. 32 the Is peak is corrected with the 0{as) corrections for the position of the pole 
and the corrections of the wave functions |jl2| (without correction to the continuum wave 
functions). Sohd hues are calculated within the summation approach , dashed lines with 
the Green function approach for nit = 180 and an electron polarisation of 0.6. The upper 
curves correspond to a pure Coulomb potential. The lower curves are from a potential 
including only the one loop correction. The dashed-dotted curve is the renormalisation 
group improved result. 

The mass shift for the levels can be obtained from the relevant part of the potential 
by performing the Fourier transformation into coordinate space, where the integrations 
can be done analytically . The surprisingly simple result is 



(3.46) 



167™^ am' 

where is the n-th logarithmic derivation of the gamma function and 7 denotes Euler's 
constant. 

The picture shows that this approach would give a good result for the position of the 
rising edge of the cross section as well as for the position of the Is peak if it were not 
necessary to sum up the leading logs. 

In view of the above mentioned complications we will always include the resummed 
potential in the following. 

Let us now return to the bound state correction for the toponium width. 
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As can be seen by inspection of eq. |2.148| we could have obtained the correct results for 



the level shifts simply by calculating the expectation value p.l48| between a wave function 
of a Schrodinger equation with a QCD potential. The term 



Hi = iV 



p' 



2m2 



(3.47) 



can thus be viewed as an absorptive part of the Hamiltonian. Adding this term is equiv- 



alent to the introduction of a "momentum dependent width" |^5|, 0. It has been shown 



45| , |51| that the 0{as) corrections to the absorptive part of the Hamiltonian vanish. 

Following this philosophy to estimate the effect of the bound state corrections to the 
width we use the Hamiltonian 



m 



2m^ 



QCD\ 



(3.48) 



where Vqcd is the one loop resummed potential. 

The results of our numerical studies for rrit = 180, p 



are shown in fig. 3.3 and 3.4 




E/GeV 



Fig. 3.3: Effect of the boundstate corrections to the toponium decay width on the total 

lold in units of a^^. Full line: without 
state corrections to the decay width. 



cross section near threshold in units of o"^^. Full line: without, dashed line with bound 



63 



R 



E/GeV 

Fig. 3.4: The same as Fig. p.3| with a different scahng - demonstrating the smallness of 

the effect. 

The cross section is shghtly enhanced and the effect is of the order of less than 1% and 
thus "really" an O(a^) effect. However, we would like to indicate that this enhancement 
is not only at the position of the IS resonance as it would be expected from a simple 
summation over resonances and from the model calculation in ref. Furthermore, it 

is interesting to note that the Hamiltonian |3.48 gives the same result as a similar one 
with the second term missing but with the replacement m —>■ m — iV /2. This convinces 
us that Hi can be considered as small. 

While the above result in the Green function approach shows reasonable results, we 
nevertheless want to add some words of caution. One observation concerns the definition 
of the perturbed resolvent. Assuming that relativistic effects are of the same order of 
magnitude and can be included in a perturbative manner, it seems necessary to go back 
to the perturbative expansion and subtract the already included parts. But due to our 
nonrelativistic approximation Hi oc p ^ the second term in the series 

ImG = ImGo + Im(Go^iG'o) + ... 

is divergent while the relativistic expression is clearly well behaved. Also the cancellation 
of gauge dependent contributions coming from the high momentum region seems to be 
more involved than it is for the level corrections. 



The second observation concerns the comparison of the results of fig |3.3| with an 
approach where one calculates corrections to the wave functions and to the position of the 
poles and afterwards takes the sum over the whole spectrum to obtain the Green function 
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46| . This, in general, cannot lead to a correct result since Hi makes the Hamiltonian non 



hermitian and thus the latter need not have a spectral representation. 

We conclude that certain electroweak and strong corrections of the same order of 
magnitude must be still missing or, even worse, it is not yet clear how to correctly add 
them. Thus a better theoretical understanding of such effects is desirable. 

However, these effects will be of the same order of magnitude and we thus may con- 
jecture that the theoretical uncertainties of the calculation of the tt cross section near 
threshold remain at the percent level. 

An effect larger than that could be induced by a relatively light Higgs boson. This 
assumption is in agreement with recent electroweak data and the corresponding fit to the 
MSSM. If the MSSM and the Higgs mechanism turn out to be realized in nature the 
lightest higgs has to lie in the range 90 to 130 GeV. Recent numerical studies show 
an enhancement of the cross section at the 10% level for itih lOOGeV which seems to 
be the value presently favoured by the analysis of radiative corrections to the SM and 
experimental limits. 
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3.1.2 Forward- Backward Asymmetry 

In the last subsection we investigated the total cross section for ti production in a future 
e'^e~ collider. This should be rather good observable to determine the mass of the top 
quark. But since the total cross section at threshold is also significantly determined by 
other quantities like the strong coupling constant and the top decay-width, it is necessary 
to have some independent observables. It was first proposed in |^ that the forward back- 
ward (FB) asymmetry would provide such a quantity. It has the additional advantage 
that it is independent of the absolute normalization of the cross section measurement 
and is therefore a quantity which can be determined experimentally with high accuracy. 
However, treating this problem in a nonrelativistic context, one encounters unphysical 
divergences which made necessary the introduction of an unnatural cut-off. In this sub- 
section we will show a systematic way how to avoid such divergencies. 

The FB asymmetry is defined as the relative difference of the cross section for particles 
produced in the forward and in the backward direction, with respect to the direction of 
the e~ beam. 



ApB ■= = (3.49) 

tot cr tot 

Since the top quarks will be identified by their main decay products W and b let 
us consider the cross section for the process e~^e~ —>■ ti bW^bW~. Production of 



W bW b via other channels can be treated as in 60 



da,+,-^hw+iw- = ^\Ma + Mv\^d<^w+bd<^w--b{'^'^yS{ti + t2 - h - k2 - bi - fealS.SO) 
with 

d<^w^b = 4\-,4^,{27r)e{ho)S{e - M'){2n)eib,o)S{bl) (3.51) 

and similarly for the decay of the antitop. 

The FB asymmetry originates in the interference terms of the vector coupling (My) 



and the axial coupling {Ma}- This is depicted in figure |375| where V indicates a vector and 
A an axial vector coupling, the star denotes complex conjugation. The matrix element 
M is essentially the vertex 'yti and Zti. Near threshold the perturbative treatment of 
the Coulomb interaction is no longer valid. Instead one has to include the whole rung of 
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Fig. 3.5 



Coulomb interactions and even worse, it is necessary to include the leading logarithmic 
contributions from the gluonic self energy corrections, as we have seen in the last section. 
We can sum the relevant set of graphs by means of the equation 

r'^ = - iKDrf" (3.52) 



for the vertex function F''. Comparison with the BS-Equation |2.26| shows that it has the 
solution 

= -tD-^GY (3.53) 
F:; = -iD-'Gj,r, (3.54) 

where it is understood, that the two legs of G on the right hand side are connected with 
the 7-matrices and a momentum integration is performed. 
We can split the matrix element into 

M\ = mim2DT\ (3.55) 

where 



mi = e^^\h)u,{h,){-i . V;,rP-) (3.56) 

V2 sm Uw 

= et'\h){-z^^^—V,sl''P-)vy{b2). (3.57) 
V 2 sm &w 

Near threshold the propagators in D and the Green function G should be dominated by 
small momenta and therefore it should be possible to replace them by their nonrelativistic 
approximations. However, for D the approximation to 0{a) is needed for the axial vertex, 
since the zero order propagators would give zero due to A~757*A"^ = 0. 
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D ^ (3.5^ 



This gives rise to the effective nonrelativistic axial vertex 

k 

757' ^^[7N7l75A-. (3.60) 

We note aheady here that the term |^ introduces an additional power in p which is not 
present in the relativistic propagator. Instead the factor ^ in |3.60| would be replaced by 

k 

which is finite for p 00. 

Using the fact that the heavy quarks are on shell up to O(a^) we can write 

j d^w+hY.\^^^ = -2ImS(m^)A+ + 0{a^) (3.61) 

(T,A 

in agreement with section 2.5. 

From the modulus squared of the propagators we obtain a factor 



dpo. 1 



(3.62) 



uj = - — [p"^ — mE — imT) (3.63) 



Collecting everything from above, performing the trace in a frame where the leptons 
move along the z-axis, and using eq. |3.16| we get for the FB-asymmetry 



o"FB,A = {cav + dAvp)^pf^ / p^dpdQFBRe[G*{p) / -——rG{p,q) ]cr^^(3.64) 

TT^F"^ Jo J [ZTT)'^ m 



where 



CAv = am^'^ + a^'^W^ (3.65) 
dAv = a(^)a(^) + (3.66) 

J dQpB = J^" d^{j^ ' dO - d6) sin (3.67) 
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The SM values for a^^\b'-^^ are given in |3.9| and p.lO| . We introduced an intermediate 



cut-off A in p.64| since the p integration is logarithmically divergent. It will be removed 
in the following. The Green function G{p) in |3.64| is defined by 



G*ip) 



(2n) 



^ Arc f°° 
:G*{p,q) = — / drr sin prGi=Q{r,0). 



p Jo 



(3.68) 



Gi=o{r, f') denotes the S-wave Green function in configuration space, eq. |3.36| . Let us now 
investigate the term 



^G{p,q)^— = ^ f d'xe-^^^^i-G{xMy^o 



(27r 



m 



m 



(3.69) 



Using the representation |3.26| one can show that only P-waves (/ = l,m = 0) contribute 
to the sum. Thus we need to evaluate eq. |3.28| for / = 1. The regular and singular 
solutions behave at the origin as 



(3.70) 



r 



respectively. We now construct a singular solution with c> = 1 out of two arbitrary 
solutions Upi and Up2 : 



ap[up2ir) - BpUpi{r)] 



by means of 

Bp is determined by the requirement of regularity at infinity 



Up := limr ^[Up2{r) - BpUpi{r)] \ 



r^^Upi{r) 



(3.71) 
(3.72) 

(3.73) 



Then the condition |3.32| demands c< = m/3. Performing the different at ion in the z- 
direction in the limit y — *• in p.69| leads to 



/ 



j—r^G{p,q) 



cos 6 



p"^ Jo 



°° smpr 

dr[ p cos pr)gy [r) 



(3.74) 



While the expressions |3.68| and |3.74| are well defined, it has been observed in [^] and 
|57| that p.64| is logarithmically divergent. Thus cutoffs (A) have been introduced to 
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make numerical predictions possible. Since in our present case the divergence is only 
logarithmical, this leads to phenomenologically reasonable results. But the explanations 
given in these references are different and so are the cut-offs. It may, however, happen that 
for the comparison with experimental data with some cuts a reintroduction of some kind 
of cut-off will be necessary 0, but we feel saver by first giving a satisfactory theoretically 
prediction and leaving this possibility open to the experimentalists. 

In any case, the introduction of these cut-offs leads to some numerical uncertainty in 
the result, and is clearly unsatisfactory from the theoretical point of view. Furthermore 
we will see in the next section, that a similar divergence occurs in the context of the 
stop-anti-stop production. Since this divergence is linear, a better understanding of the 
origin of this kind of divergencies is necessary to give a quantitative prediction. 

The key point of our solution of the problem is to go back to the perturbative sum of 



the graphs for (Tfb- This is illustrated in figure |3T5 . 

Consider now the leading (tree) contribution in the nonrelativistic approximation (e.g. 
the first graph in fig. 3^). In this case both S- and P-wave Green functions are replaced 

GiP,p-) ^ <|M^. (3.T5) 



m 



This leads to the logarithmically divergent expression (for A ^ cxo) 



^^^'^ = J ^^^^ Jo ^^ |£-E-.rp - ^^-^^^ 



m 



Since the remaining graphs in fig. by power counting can be shown to give finite results 
we conclude that the nonrelativistic approximation was not valid in the tree graph due to 
the extra power in p from the axial vertex. Furthermore since 3.76 is the only divergent 
contribution in the (infinite, but assumed to be convergent) sum |3]^ representing G we 
conclude that the divergent part of |3.64| is entirely contained in the first (free) contribution 
3.76| . We could now remove this divergence by a replacement p/m ^ p/Ep, as indicated 
above. This, however, does not respect another qualitative difference between the exact 
tree contribution and its nonrelativistic approximation. Namely in the relativistic calcu- 
lation the phase space is cut off when the momentum squared of one quark falls below 
the invariant mass of the decay products . 



Therefore, we return, instead, to the relativistic expression for the tree contribution. 
Since we are dealing with a decaying particle, we will have to include the self energy 
^Thc author is grateful to M.Jezabek for this information 
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contribution due to the decay in the propagator. This will in general lead to gauge 
dependent results, but we can use the constant on-shell width to obtain the leading 
contribution in the weak coupling 

This makes the quantities we are considering finite (for F ^ 0) and gauge independent 
and the remaining higher order contributions (e.g. gauge dependent if one considers 
only one diagram) can be calculated perturbatively. Therefore to leading order it seems 
reasonable to take the constant width approximation for the top propagator. 

Since we are only investigating processes with ti intermediate states which will give 



the main contribution to the cross sections we focus on the graph shown in fig. 




Fig. 3.6 

A further advantage of this method is that now the other - non resonant - graphs for 
the process e'^e~ W^bW~b to leading order in the weak coupling need only to be 
added to yield the background contribution. 

Performing the straightforward calculation for the relativistic tree contribution with 
a constant, but non-zero width F, we arrive at 



cr 



(1) 

FB 



p2 



where 



A(p^ 



mF 



This could have been also obtained by replacing the fermion propagator by 

F (;^ + /i) 




S 



djj 



7r[(m2 — /i^)2 + m^F^] — fi^ + ie 
and cutting it to effectively replace the phase space element for a stable quark 



(3.79) 



d<!> 



stable 



(2vr) 



;Q{po)S{p -m 



71 



by 

d<^unstable = J^^Q {po) Hv") ■ (3.80) 

To conclude we can say that our simple prescription to obtain finite, gauge indepen- 
dent results to the desired order is to replace the nonrelativistic tree contribution by the 
relativistic one and leave the (finite) rest unchanged: 

(^FB = (^FB + {crFB,A " (^FbT) (3-81) 

A^oo 



In fig. p.7| the numerical results are shown for the different contributions to apB for 
A = 300 GeV. This can be estimated to give a result for limA^oo(o"F_B,A — (^fb^a) lyi^g 
only 2% below the final answer. The results of the purely nonrelativistic calculations are 



compared to our approach in Fig. |3.8| for nit = 180 and an electron polarization of 0.6. 



We also included the hard corrections as given in |57 |. 

We believe that our approach has the decisive advantage that it can be based upon to 
the original set of graphs to be considered and it is clear which graph has been calculated 
to which accuracy. Thus at least in principle a systematic improvement is possible. The 
present approach should also be applicable to the 0{as) final state corrections calculated 
in |2|. 
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Fig. 3.7: Different Contributions to apB 
dashed: cr|?g, dotted:cr^^^"''', dashed-dotted: ap^^ 
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3.2 Axial Contribution to the tt Total Cross Sec- 
tion and the Production of Stop-Antistop near 
Threshold 

The total cross section for the production of P-wave states near threshold gives the lead- 
ing term for the production of stop-antistop and a contribution of 0(0;^) to the total 
cross section for ti. Only a qualitative investigation on the basis of the Coulomb Green 
function for the former has been published up to now |61]. We will apply the method 



of replacing divergent nonrelativistic graphs, implicit in the Green function approach, by 
the relativistic, finite ones to get a quantitative reliable result. 

Considering the Coulomb green function as a first approximation it was observed 
in 1^ that the total cross section for a stop-antistop pair near threshold develops an 



unphysical linear divergence due to the nonrelativistic approximation. Even a relativistic 
calculation of the imaginary part of the one loop contribution using the commonly used 
propagator — -|- irriT) remains logarithmically divergent. This is due to the fact 

that this propagator does not fulfill the requirements of quantum field theory as does the 
propagator 

dfi^ niT 1 ^ 

which is in agreement with the Lehmann representation of the full propagator. The tree 
contribution using this propagator reads 



whereas we have for the nonrelativistic tree contribution 

3r , , . /-^ , 



p2 ' p4 

(3.83) 



a 



1,"'" 



+ dsp) / dp— ^o-^/i- (3-84) 

However, a single subtraction of the tree graph is not sufficient. One should also subtract 
the logarithmically divergent one-gluon exchange term. But since we replaced the gluon 
propagator by the resummed one we should also calculate the relativistic graph with a 
resummed propagator to get the correct behavior at infinity. Unfortunately this leads 
to renormalon ambiguities. Therefore, we conclude that the uncertainties in choosing 
the right cut-off are of the same order of magnitude as higher QCD corrections, and we 
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circumvent these difficulties for the time being by keeping the cut-off in the logarithmic 
divergent terms, choosing its value A = Aq = m from the observation that this cut-off 
would have given the correct answer in the case of the FB-asymmetry discussed in the 
last section. 



a 



0.04 



0.03 



0J)2 



0.01 



E/GeV 



Fig. 3.9: Different Contributions to 
dashed:(j^°,dotted:crf5^r) dashed-dotted: ai^rei 



The same is true for the pure axial contribution to the total cross section for tt production 
since each axial vertex contributes in the nonrelativistic limit an extra power in p as 
explained in the last section. The only change is to replace the factor 

3 



(3.85) 



by 



ff:=6{cA + dAp). (3.86) 
Using the methods of the preceding sections we get within the Green function approach 



1 



3T 



X 



dpp^ / dn 



d^q q^^^ 
(27r)3 m 



G{q,p) 



a 



fifii 



(3.87) 



where fx denotes either fs or The different contributions to the cross section for the 
scalar case 

a = at°-at°,, + a^ (3.88) 
are depicted in fig. ^]9| for cos^ 9t = 0.5, m = 180, F = F^p. The net result is shown in fig 



3.10| for the axial contribution to the tt cross section. 
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R 
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E/GeV 
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Fig. 3.10: Drell-ratio (R= (Ta/ct^p) for the production of tt 



According to ref. ||62| we included the hard corrections by a factor (1 — ^)^. We 
conclude that the axial contribution gives a sizeable effect (of a few percent ) to atot 
beginning at ^ 5GeV above threshold. However, it is small enough in order to hide the 
uncertainty in A with respect to O(a^) corrections to the vector contribution to atot- 



Finally in fig. |3.11| we compare the cross sections for the production of tt near threshold 



for two different values of the decay constant ( cos^ 6t = 0.5, m = 180). 



0.03 




E/GeV 



-3-2-10123 

Fig. 3.11: dashed curve: F = l.SGeV, full line: F = O.SGeV 
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Chapter 4 
Conclusion 



In this thesis we considered the quantum field theory near thresholds of ultra-heavy par- 
ticles, where we especially focused on the tt and it systems. Due to the large mass and 
the large width of the top quark it is possible to treat the former in a pure perturbative 
manner although the top is a strongly interacting particle. The large mass of the top 
close to the weak scale may be an indicator of new physics. Therefore, it is of paramount 
interest to predict physical observables within the framework of the existing Standard 
Model in order to be able to single out possible new effects. 

In this context the question of the possible large corrections due to a running width of 
the top quark could be solved in this work. It was shown that a previously not considered 
contribution leads to large cancellations and leads especially to a gauge independent result. 
This was first achieved in the narrow width approximation and subsequently improved by 
the introduction of a new relativistic zero order equation already including the on-shell 
decay width. With the help of that equation and a Ward identity it was further possible 
to show that the bound state corrections to the decay width for a generic fermionic system 
which does not decay by annihilation is given by 



Furthermore a systematic derivation of the tt potential to numerical order 0(q;^) was 
given. Especially the contribution of a QCD box graph and a thorough investigation of 
effects of the weak interaction are new. Also the potential for scalar-scalar bound states 
has been investigated with the help of a new zero order equation for such systems. A 
deviation from a previous result was found. 
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In the second part of this work we considered within the numerical Green function 
approach the cross sections for the production of ultraheavy particles near threshold. 
First we derive the general formalism and consider the total cross section. Here the 
large width is demonstrated explicitly to be very important for hiding nonperturbative 
effects. Then we investigate the forward-backward asymmetry near threshold. We avoid 
the introduction of an unnatural cut-off as contained in previous work. This is achieved 
by identifying the divergent graphs in the nonrelativistic Green function and replacing 
them by the relativistic ones. The same method can also be applied to the pure axial 
contribution to the ti total cross section as well as to the stop-antistop production. This 
enables us to give first quantitative predictions of these observables. 

While some aspects of the production of heavy particles near threshold are well under- 
stood by now, there remain several unsolved problems. Especially proceeding to higher 
orders in perturbation theory seems rather difficult in view of the necessity of including 
the leading Coulomb interaction to all orders. It has not been possible to completely 
combine the rigorous bound state approach of Chapter 2 with the practically tractable 
numerical approach of Chapter 3 in a satisfactory way, yet. Nevertheless, we hope that 
this work provides a basis for further investigations in this direction. 

ae 
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Appendix A 
Expectation Values 



In sect. 2.4 and 3.1 we needed the expectation values of logarithmic potentials between 
Schrodinger wave functions. They can be obtained by 



,ln"r 



dX 



(A.l) 



if the expectation value {r^ ^) is known analytically. In the latter the representation 



Jn— 1 



- / - 1)! c/2"-'-i 



:i-z) 



(A.2) 



of the Laguerre polynomials may be used. This allows an easy evaluation of the integra- 
tions and the remaining differentiations can be done with some care afterwards: 



n-l-1 / _\ ^ 

,n — I — 1 — k J 



k=0 



-21-2- 
k 



Using now the Fourier transformations 



In^ 

q ^ 
In^a; 

q ^ 



7 + In /ir 
27rr 

1 TT^ 

■[— + 2(7 + ln/ir)2] 



27rr 6 



(A.3) 

(A.4) 
(A.5) 



one arrives immediately at eq. p.46| and |2.83| , respectively. 

38 
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